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ZENER  OSCILLATIONS 


1.  INTRODUCTION 

The  question  of  the  existence  of  Zener  oscillations  has  remained 
controversial  for  more  than  forty  years.  Experimental  observations 
indirectly  supporting  their  existence  have  been  reported  by  Koss  and 
Lambert  (1)*,  but  so  far  this  widely  quoted  work  is  the  only  credible 
experimental  evidence. 

In  the  present  report  we  continue  the  Zener  oscillation  studies 
begun  in  the  1980  progress  report  on  Semiconductor  Millimeter  Wave¬ 
length  Electronics  (2) .  In  that  report  an  extensive  critical  review 
of  band  structure  dominated  carrier  dynamics  was  presented.  It  was 
concluded  that  the  phenomena  limiting  the  realization  of  Zener  oscilla¬ 
tions  are  scattering  and  Interband  tunneling.  A  discussion  and  numeri¬ 
cal  estimate  of  the  tunneling  probability  indicated  that  this  is  not  a 
serious  limiting  factor,  and  it  was  concluded  that  Zener  oscillations 
would  require  an  adequate  lengthening  of  the  scattering  lifetime. 

The  motivation  for  studying  Zener  oscillations  is  their  potential 
of  realizing  a  device  for  generating  tunable  submillimeter  radiation. 
Like  any  other  effective  electronic  generator  of  radiation,  a  practical 

*The  numbers  in  parentheses  in  the  text  indicate  references  in  the 
Bibliography. 
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Zener  oscillator  will  have  to  rely  on  phase  coherence  of  the  electrons 
coupled  to  the  radiation.  A  possible  scheme  for  obtaining  the  required 
phase  initialization  and  phase  focusing  has  been  suggested  by  D.  L. 

Rode  (private  communication)  and  will  be  reported  on  elsewhere.  In 
the  present  work  we  discuss  the  theoretical  aspects  of  such  a  device 
from  the  viewpoint  of  quantum  theory.  Specifically  we  examine  an  elec¬ 
tronic  state  analogous  to  the  coherent  states  of  a  harmonic  oscillator 
("Glauber  states" (3))  with  a  wave  function  whose  mean  square  position 
and  momentum  uncertainty  product  approaches  the  minimum  uncertainty 
level.  Such  a  state  can  be  described  by  a  wave  packet  of  states  of 
the  conventional  representation.  The  dynamics  of  such  a  wave  packet 
must  be  determined. 

In  view  of  our  goals  we  have  inquired  into  the  physical  meaning  of 
several  existing  solutions  of  the  problem  of  an  electron  in  a  crystal 
in  the  presence  of  an  applied  electric  field.  We  begin  by  quoting  the 
standard  theoretical  methods  and  their  results.  No  derivations  are 
given,  since  they  can  be  found  in  the  1980  report  and  in  references 
(A),  (5),  (6)  and  (7). 

Based  on  the  interpretation  of  these  results,  we  propose  two  kinds 
of  wave  packets  to  represent  coherent  band -electron  states.  Further¬ 
more  we  calculate  absorption  and  emission  probabilities  of  Zener  os¬ 
cillations  in  a  one-band  scheme.  These  probabilities  are  directly 
related  to  the  spectral  analysis  of  the  Zener  oscillations  in  a  given 
band  structure.  We  report  on  a  computation  of  these  spectral  compo¬ 
nents  for  a  variety  of  applied  field  directions  for  the  conduction  band 


of  GaAs. 
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2.  THEORETICAL  BACKGROUND 


The  standard  technique  for  studying  the  electronic  properties  of 
crystals  is  the  one-electron  model.  In  such  a  model  the  electron  is 
viewed  as  moving  in  an  average  periodic  potential.  It  is  assumed  that 
in  an  applied  electric  field  this  average  potential  is  unaffected. 

This  approximation  is  based  on  the  relative  weakness  of  the  applied 
field  when  compared  with  the  effective  fields  of  the  lattice  bonds. 

There  are  several  different  approaches  to  the  solution  of  the  one 
electron  model  equations  in  the  presence  of  an  applied  electric  field. 
In  all  cases  the  Schroedinger  equation  to  be  solved  is 

,t)  (1) 

(choosing  the  sign  convention  e  *  -|e|). 

For  our  purpose  it  is  useful  to  show  the  solution  in  each  of  three 
different  representations:  The  Bloch,  or  crystal  momentum,  Wannier  or 
lattice  site;  and  the  kq  representation. 

2.1  FORMULATION  OF  THE  SCHROEDINGER  EQUATION  IN  THREE  REPRESENTATIONS 

2.1.1  Bloch  Functions 

A  natural  basis  for  crystal  electron  wave  functions  is  the  set  of 

Bloch  functions  \p  (k,r).  The  Bloch  functions  are  eigenfunctions  of  the 
n  -2 

periodic  Hamiltonian  Hq  M  ^  +  V(r),  and  they  are  labeled  by  the  two 
Indices  k,n.  The  index  n  is  the  band  index  defined  by  the  Hamiltonian 

Vn(^)  -  En(k)<|>n(k,r) . 

The  index  k  is  defined  by  the  lattice  translation  operator  T(am)  * 
exp (ip *a ^  through 


(£. 
\  2m 


») 


+  V(r)-  ee*r  j  4»(r,t)  =  i  -ft 


T(am)^n(k,r)  -  ei^k‘am^n(k,r). 

The  Indices  n,  k  are  good  quantum  numbers  since  they  are  generated  by 
commuting  operators 

tTUm),  H0]  -  0. 

The  presence  of  the  field  term  ee*r  in  Equation  (1)  breaks  the  period! 
city  of  the  Hamiltonian  and  k  is  no  longer  a  good  quantum  number. 
Therefore  we  write  the  solution  to  Equation  (1)  in  the  Bloch  function 
representation  as 

i|>(r,t)  -£/  dk<}>n (k , t ) (k,r)  . 
n 

If  we  choose  the  x  axis  to  lie  along  the  field  e,  substitute  in  the 
Schroedinger  Equation  (1) ,  and  evaluate  the  matrix  element  of  the  posi 
tion  x  that  appears  in  the  field  term  (Appendix  6.3) 

<n,k|x|n' ,k’>  =  i5nn,  g|—  6<k-k')  +  X^.SCk-k')  (2) 

x 

we  obtain  the  equation  of  motion  of  the  envelope 

IEn®  - lec  w  - 1 11  k  iv*'c)  -  «£  y/'1’  <3> 

x  e 

whece  Xne  '  IT  %S’?>d7 

is  the  polarization  matrix.  The  integral  is  taken  over  a  unit  cell, 
with  n  the  volume  of  the  cell. 

2.1.2  Wannier  Functions 

Another  useful  representation  is  formed  from  the  set  of  Wannier 
functions,  defined  in  terms  of  the  Bloch  functions  by 
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W  (r-R  )  =  /  °  -  f 

n  n  V  <2*)3  B.Z. 


— ik*R - — 

e  n^n(k,r)dk 


where  R  is  a  lattice  site  vectr r. 
n 

In  the  Wannier  representation  the  solufon  of  Equation  (1)  is 


<Kr,t)  -EE  b  <R  ,t)W  (r-R  ) 
n  m  n  m 

n  m 


where  b  (R  ,t)  is  a  solution  of 
n  m 


[E (-iV)-i  h  3-  2  _  b  <r,t)  -EE  Unn'(Rm'»R  )bn'(Rm 
n  dt  _  —  n  ~  •=*  nn  m  m  n  m 

r  Rm  n*  R  , 

m  m 


,,t)  = 


0 

(A) 


where  U  , (R, ,R)  -  f  W  , (r-R  , )ee*rW  (r-R  )dr 
2.1.3  The  kq  Representation 

A  representation  specifically  designed  for  a  periodic  medium  which 

suppresses  the  band  index  (its  basis  functions,  in  terms  of  the  Bloch 

or  Wannier  representations,  are  sums  over  all  the  bands)  is  generated 

by  the  translation  operators  in  the  direct  and  reciprocal  lattice. 

These  operators  T(a  )  -  exp (ip ’a  )  and  T(b  )  ■  exp(iq*b  )  are  a  com- 

n  n  m  m 

plete  set  of  observables  [a  ,b  ]  ■  2T5  and  they  define  wave  functions 

n  m  nm 


I(V*M  ’  exp 


T<W  ■ eitp 


(k  ^*an)^  kq^ 

(*«*.)»* 


with 


^  'J  (2,)  3  T 


exp (ik'R^) 6 (r-q-R^) 


where  a  is  a  lattice  vector. 
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In  the  kq  representation  the  solution  of  Equation  (1)  is 

i|/(r,t)  =  f  dk  f  dqC(k,q,t)^-(r) 

where  C(k,q,t)  is  a  solution  of 

^  +  VC«1>  —  ce'-Ci  -l-q)JcOc,q,t>  =  ±ii  |^-  CJClc.q.t)  (5). 

2.2  THE  ACCELERATION  THEOREM  IN  THREE  REPRESENTATIONS 

The  solution  of  the  Schroedinger  equation  can  be  carried  out  in 
each  of  these  representations.  The  full  solutions  must,  of  course,  be 
identical,  but  each  affords  a  different  view  of  the  problem.  This  will 
enable  us  to  gain  physical  insight,  to  choose  superpositions  of  states 
to  correspond  to  various  initial  and  boundary  conditions,  and  to  make 
appropriate  approximations. 

The  solution  of  Equation  (1)  and  of  its  equivalent  forms  Equations 
(3),  (4)  and  (5)  is  the  motion  of  an  electron  in  an  electric  field, 
thus  an  acceleration  theorem.  We  shall  take  up  the  solution  in  the 
three  representations  in  the  reverse  order  to  their  introduction, 
because  in  that  way  we  can  proceed  from  the  most  general  and  abstract 
to  the  most  intuitive  form  of  the  acceleration  theorem. 

2.2.1  The  kq  Representation 

In  this  representation,  generated  by  the  translation  operators  in 
direct  and  reciprocal  space,  the  coordinate  and  momentum  operators  are 


These  operators  are  defined  within  the  unit  cell  of  their  respec¬ 
tive  spaces,  but  because  the  translation  operators  contain  an  arbitrary 
phase  factor  2niT,  they  are  not  localized  in  a  particular  cell  and  act 


equally  in  all  the  cells  of  their  space.  Another  way  of  stating  this 


property  is  to  say  that  fc-—  places  the  electron  on  an  infinite  point 
lattice  in  both  direct  and  reciprocal  space,  with  the  lattice  point 
exactly  localized  both  in  the  unit  cell  and  in  the  Brillouin  zone, 
does  not  belong  to  a  band;  indeed  the  exact  localization  within  the  cell 
is  possible,  in  the  language  of  the  more  familiar  Bloch  or  Wannier 
representations,  because  contains  a  superposition  of  all  the  bands. 

For  this  same  reason,  the  acceleration  theorem,  deduced  by  Zak  in 
the  form  of  a  Heisenberg  operator  equation  of  motion 


f -|[H,k]  -^ee 


is  exact  and  not  a  one-band  approximation.  It  is,  of  course,  this 
property  that  prompted  Zak  to  invent  this  representation. 


The  Heisenberg  operator  k(t) 

;  S  H<t-C0>-  -  5 

k  *  e  ke 


is 


The  physically  meaningful  crystal  momentum  is  the  expectation  value 
«Kr,t0)  |k(t)  |iKr,tg)>  *  <ip (r , t)  | k | (r , t) >  . 

2.2.2  Wannier  Representation 

The  Wannier  representation  leads  very  directly  to  a  very  useful 
correspondence  between  classical  and  quantum  dynamics  of  the  crystal 
electron  quasiparticle.  In  this  representation  k  is  an  operator 
(i  grad  — ) •  If  the  interband  and  nonlocal  terms  n  /  n' ,  m  4  m'  are 
dropped  in  Equation  (3) ,  the  one-band  Schroedinger  equation  for  "weak" 


fields  is  reduced  to 
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[E  (-iV)-i  ft  +  ee*r]b  (r,t)  =  0  . 

n  ot  n 

Here  E^C-iV)  is  the  modified  kinetic  energy  operator  that  includes  the 
effect  of  the  periodic  crystal  potential,  and  the  envelope  function 
bQ(r,t)  can  be  viewed  as  describing  a  wave  packet  of  Wannier  functions. 
Its  trajectory  is  given  by  the  equivalent  classical  Hamiltonian 


E  (-iV)  +  ee‘r  E  (k)  +  ee*r 
n  n 


This  correspondence  recognizes  ft  k  as  a  crystal  momentum,  and  leads  to 
the  classical  Hamilton  equations  of  motion 

^  3H  13V» 

r  3p  ft  3lc 


—  1  Oil 

p  -  ft  k  -  -  -r= 


3H  - 
"  3r  “  e£ 


which  can  be  integrated  to  give  the  classical  position  of  the  electron 


r (t>  -r(tQ)  *  “  lEn(k(t))  -En(k(t0))] 


(6) 


with  k(t)  -  k(tQ)  -  (t-tQ). 

It  should  be  remembered  that  this  elegant  deduction  of  the  Zener 
oscillation  dynamics  is  based  on  a  wavepacket  formalism  in  a  one  band 
scheme . 

2.2.3  The  Bloch  Representation 

The  Bloch,  or  momentum  representation  leads  very  directly  to  a 
group  of  intuitively  appealing  results  on  Zener  oscillations.  Starting 
with  Equation  (3)  we  can  obtain  an  exact  equation  of  motion  for  the 
probability  density  of  the  momentum  distribution 
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(ee  dr + *  &  2>n°° 3  0 

x  n 

which  is  satisfied  by  any  arbitrary  initial  distribution  G(k)  that 
changes  with  time  according  to 

ZUn(k)|2  -  G(kx-ee  ky,  kz). 
n 

Hence  for  any  initial  superposition  of  momentum  eigenstates,  the  expec¬ 
tation  value  of  the  momentum,  defined  as 

<k>  «£  /  |*n(k)|2k  dk 
n 

behaves  as 

<kx(t)>»<kx(t0)>+f  (t-t0)  • 

Note  again,  as  in  the  kq  representation,  that  this  acceleration  theorem 
is  obtained  without  neglecting  interband  terms;  unlike  the  classical 
trajectory  of  Equation  (6),  the  time  dependence  of  the  crystal  momen¬ 
tum  is  unaffected  by  interband  mixing. 

In  a  one-band  model  Equation  (3)  becomes 

>E„®  -  etXnn-lee  IT  *  1  "  IT  • 

We  can  write,  for  the  polarized  band  structure 
E^X)(k)  -  E(k)  -  eeXnn; 

Then  the  solution  can  be  written 

_  -  i  E  (t) 

♦n(k,t)  -  (k)e  ”  v 


ft" 


1 


ft 


•  < 


'  ■  I 

ft  ' 


ft 


with 


-10- 


Vk)  = j=  6  Vkyo)6(kz-kzo)exp  [if/  X(VE(1>(k))dkx]  (7) 

where  <  is  the  length  of  the  reciprocal  lattice  vector  lying  along  the 
field  direction  (chosen  as  the  x  axis).  (If  the  field  direction  is  a 
principal  lattice  direction,  <  is  the  "diameter"  of  the  Brillouin  zone 
along  the  x-axis). 

The  wave  function  4>v  is  periodic  in  the  extended  Brillouin  zone, 
traversing  a  phase  shift  of  2irv  radians  as  k^  traverses  <  wave  numbers. 
Thus  the  energy  eigenvalues  are 

E  -  2ttv  —  +  -  f  E(1)(k)dk  . 
v  <  <  JQ  x 

These  energies  have  been  called  Wannier  levels  or  Stark  levels, 

and  they  form  the  "Stark  ladder"  which  has  been  the  subject  of  so  much 

controversy.  We  use  the  terminology  Stark  levels  for  Ey  and  Stark 

functions  for  ,  to  avoid  confusion  with  the  Wannier  functions 

W(7-R  ). 
n 

It  is  possible  to  form  a  superposition  of  all  the  Stark  states  of 
a  crystal  with  equal  weight  for  each  state  in  such  a  way  that  the  re¬ 
sulting  superposition  has  a  def inite-albeit  time-dependent-value  of  kx» 
We  only  indicate  the  result  without  giving  any  of  the  intermediate 
manipulations : 

>Kr»t)  /  4>v(k)e  V  <Kk,r)dk 

V  ~ 

est 

■  Et,kyo,k20)  *  dkxE(1)(k)j  (8) 

where 
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*(l  “•  Vk2o)  '  ua0t(t)  >r>el  k<t)'r 

is  a  "time  dependent  Bloch  function".  The  designation  "Houston  func¬ 
tion"  is  used  for  both  4>(r,t)  and  iKk(t),r). 

2.2.4  Summary 

By  reviewing  these  three  approaches  to  the  acceleration  theorem  we 
have  attempted  to  emphasize  the  relation  of  the  dynamical  and  geometri¬ 
cal  aspects  of  the  electron  motion  in  a  crystal.  The  time  dependence 
of  the  crystal  momentum  is  an  exact  result,  independent  of  the  one-band 
approximation.  The  Stark  ladder  and  the  localization  of  the  electron 
in  a  correspondence  principle  sense  does  require  this  approximation. 

This  must  be  distinguished  from  the  localization  in  the  kq  representa¬ 
tion  which  has  no  classical  analogy.  We  will  study  questions  of  charge 
localization  below,  after  a  brief  mention  of  the  current  literature. 

2.3  RECENT  LITERATURE 

The  quaslclassical  electron  dynamics,  Zener  oscillations,  and  their 

quantization  in  a  Stark  ladder  of  energy  levels  are  seen  to  follow  in 

straightforward  fashion  from  a  one-band  scheme.  The  controversy  about 

the  observability  of  Zener  oscillations  or  the  Stark  ladder  revolves 

about  the  justifiability  of  neglecting  the  interband  terms  X  ,  . 

m,  nr  m 

The  Stark  ladder  could  be  destroyed  by  the  broadening  of  the  Stark 
levels  due  to  the  finite  lifetime  of  the  states. 

This  question  has  been  debated  in  the  literature  for  more  than 
forty  years.  The  only  credible  observation  that  has  been  reported  (1) 
is  the  existence  of  a  "staircase"  modulation  of  the  Franz-Keldysh  ef¬ 
fect  as  a  function  of  the  field.  This  modulation  effect  of  the  field 
on  the  interband  absorption  was  predicted  by  Callaway  (8),  and  while 


ft 


ft 


■4 


ft 


ft  ' 


ft 


ft 
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its  verification  of  the  Stark  ladder  is  indirect,  it  appears  difficult 
to  account  for  on  any  other  basis. 

Most  of  the  calculations  in  the  current  literature  appear  to  be  in 
agreement  that  in  crystals  with  moderate  or  wide  bandgaps  the  contri¬ 
bution  to  the  Stark  level  lifetime  due  to  the  applied  electric 
field-that  is  to  say,  the  tunneling  probability-is  small.  In  one  of 
the  most  recent  such  calculations  (7)  interband  tunneling  was  studied 
using  a  time  evolution  operator.  The  wave  function  used  in  the  calcula¬ 
tion  includes  the  tunneling  process  and  therefore  the  broadening  of  the 
Stark  levels  is  directly  computed.  Results  were  obtained  for 
nearly-free-electron  approximation  in  a  two-band  scheme  representative 
of  GaAs.  The  broadening,  calculated  to  second  order  in  the  field,  is 
found  to  be  neglibible  for  fields  up  to  10^V  cm 

The  theoretical  evidence  thus  supports  the  reality  of  the  Stark 
ladder.  It  suggests  that  the  limitation  of  the  observability  of  Zener 
oscillations  is  set  not  by  the  intrinsic  lifetimes  of  the  energy  levels 
but  is  to  be  ascribed  to  scattering.  The  one-band  approximation  may 
be  accepted  as  well  supported  by  the  current  best  estimates,  and  our 
further  discussion  will  be  based  on  it. 
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3.  PHYSICAL  INTERPRETATION:  LOCALIZATION  AND  CORRESPONDENCE 

PRINCIPLE 

3.1  THE  HARMONIC  OSCILLATOR  ANALOGY 

In  the  "moderate"  fields  under  discussion  (e<10^V  cm  ,  the 

classical  excursion  Ax  of  the  electron  oscillation  orbit  for  Zener 

oscillations  is  AE  /e,  where  AE  is  the  width  in  energy  of  the  nth 
n  n 

band.  Since  AE  is  a  few  eV,  Ax>10  ^  cm,  that  is  to  say,  many  times  the 

lattice  parameter  a.  Inspection  of  the  form  of  the  Stark  function 

^  of  Equation  (8)  shows  that  the  total  phase  shift  experienced  by  4>v 

as  k  traverses  the  Brillouin  Zone  is  of  order  Ax/a,  so  d>  is  many 
x  v 

electron  wavelengths  long. 

Thus  a  Stark  level  in  a  moderate  field  is  similar,  in  this  respect, 
to  a  coulombic  or  harmonic  oscillator  energy  eigenstate  with  a  large 
quantum  number.  A  particle  in  such  a  high  energy  eigenstate  is  local¬ 
ized  to  the  extent  of  having  an  appreciable  probability  of  being  found 
only  where  its  kinetic  energy  is  positive  (near  the  nucleus  or  near 
the  potential  minimum  respectively).  Within  this  range,  its  probability 
density  is  time-independent. 

However,  a  particle  with  this  much  average  energy  can  be  localized 
more  closely  by  forming  a  coherent  superposition  of  several  adjacent 
energy  states  in  such  a  way  that  at  some  particular  time  their  wave 
functions  all  add  in  phase  at  some  particular  point  along  the  orbits, 
and  cancel  elsewhere.  Such  a  localization  would  not  in  general  be 
expected  to  persist,  but  in  some  form  it  underlies  the  correspondence 
principle  for  the  formulation  of  classical  orbits  from  quantum  theory. 
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For  the  example  of  a  harmonic  oscillator  the  theory  of  such  coher¬ 
ent  superpositions  of  energy  eigenstates  is  highly  developed,  since  it 
can  serve  as  the  basis  of  the  quantum  coherence  theory  of  light.  Har¬ 
monic  oscillator  eigenstates  |a>  of  the  destruction  operator  a  have 
wave  functions  of  Gaussian  shape  with  a  localization  as  narrow  as  that 
of  the  oscillator  ground  state,  and  they  oscillate  in  the  quadratic 
potential  with  the  oscillator  frequency,  and  without  spreading.  Any 
narrower  localization  requires  a  wider  range  of  energy  eigenstates  |n> 
and  will  cause  the  coherent  state  to  spread  with  time.  The  states  |a> 
have  minimum  uncertainty  products  for  simultaneous  measurement  of 
position  and  momentum. 

Our  purpose  here  is  to  Initiate  a  similar  study  of  coherent  states 
of  a  crystal  electron  in  an  applied  field,  to  serve  as  models  for  the 
quasi-classical  electron  executing  Zener  oscillations.  Relying  on  the 
conclusions  outlined  in  Section  II  above,  we  will  confine  our  attention 
to  one-band  states.  We  have  available  two  types  of  one-band  wave  func¬ 
tions:  Houston  functions  and  Stark  functions,  eigenfunctions  respec¬ 
tively  of  momentum  and  of  energy.  We  work  out  the  relations  between 
the  two  types  of  wave  function,  and  the  localization  of  electrons  in 
each.  We  then  discuss  the  localization  of  electrons  in  wave  packets 
formed  from  superpositions  of  eigenstates,  and  an  approach  to  the  con¬ 
struction  of  minimum  uncertainty  wave  packets.  We  also  compute  transi¬ 
tion  probabilities  between  Stark  states,  since  they  are  related  to  the 
emission  and  absorption  of  Zener  radiation. 

3.2  HOUSTON  FUNCTIONS 

A  crystal  electron  generated  by  thermal  excitation  across  a  band- 
gap  (e.g.  phonon  absorption)  is  unlocalized,  and  so  it  is  not 
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unreasonable ,  in  the  presence  of  a  field,  to  represent  it  by  a  Houston 
function  (see  II  2c  above) 


-±  r 


E(1)(k(0)  +  )dt* 


*(r,t)  -  *(k(0)  +  ,  r)e 


fit  66  A 


*(k(0)  +  ^  ,  r)e 


kx<0)+^  E(l>(k.)dk. 

k  (0)  x 

x 


where  we  have  used  the  acceleration  theorem  tik  ■  ee;  this  explicitly 
represents  the  Houston  function  as  a  Bloch  function  with  time-dependent 
label  k  and  phase,  and  shows  that  it  can  be  labeled  by  the  initial 
crystal  momentum  k^CO). 

Although  the  electron  is  unlocalized,  in  some  sense  it  must  be, 
in  the  presence  of  the  field,  in  accelerated  motion,  and  the  Houston 
function  should  be  capable  of  describing  this  motion.  The  most  direct 
approach  to  a  description  of  the  motion,  the  computation  of  the  expec¬ 
tation  value  <x>  of  the  position,  fails  because  its  matrix  elements  are 
singular,  as  seen  from  Equation  (2).  The  divergence  arises  from  the 
S-  function  normalization  of  the  Houston  functions  which  we  can  write 
in  the  form 


<k^(0),  r|kx(0),  r>  -  6(kx(0)-k^(0)); 

to  obtain  physically  meaningful  results,  it  is  useful  to  avoid  such 
singularities  except  as  limits  of  finite  procedures. 

6-  function  normalization  is  commonly  used  for  wave  functions  of 
unlimited  extent  such  as  the  plane  waves  representing  free  particles, 
or  Bloch  functions,  or  Houston  functions. 


For  all  these  functions  one  could  argue  that  a  normalization  in  a 
finite  volume  V  permits  easy  physical  interpretation,  through  the  quan¬ 
tization  of  the  label  k  introduced  by  periodic  boundary  conditions. 

This  leads  to  a  normalization 

<k’lk>  "  6kk” 

and  expectation  values  such  as  <k' |xn[k>  are  easily  computed.  However, 
such  "box  normalization"  hampers  the  description  of  the  time-dependence 
of  k  for  accelerated  particles,  and  indeed  the  postulation  of  periodic 
boundary  conditions  is  questionable  at  best  in  the  presence  of  a  field 
which  destroys  the  postulated  equivalence  of  the  boundary  points. 
Therefore  it  is  necessary  that  we  deal  with  the  finitely 
non-no rmaliz able  infinite-crystal  Houston  functions. 

The  method  for  avoiding  divergent  expectation  values  can  be 
developed  by  analogy  with  free-particle  plane  waves.  Here  we  have 

i  I  nli  ^  Cm  a  ~ikrx  /  .  3  \  n  ikx  .n  3n  xt.  , 

<k  |x  | k>  - J  dxe  ^-i  e  =  i  — -  6(k-k') 

which  is  meaningless  for  k  ■  k' .  A  meaningful  expression  can  be  ob¬ 
tained,  however,  by  using  the  states  [k>  as  a  basis  for  the  construc¬ 
tion  of  normalizable  wave  packets. 

A  wave  packet 

<Kx,t)  -  /  f(k)<Kk)dk 


has  the  norm 


<«(x,t)|^(x,t)>  -  J  dx J  dk,f*(k')^*(k’) J  dkf(k)^(k) 
-  J  dk'/ dkf*(k')f(k)6(k-k') 


dkf*(k)f (k) 


and  so  can  be  normalized  by  normalizing  f(k);  evidently  the  same  proce¬ 
dure  will  yield  convergent  expectation  values  <^(x,t) |xn|^(x,t)>. 

We  shall  make  use  of  the  Gaussian  distribution 


Vk-V 


(k-k0r 

4  a2 


normalized  to 

f  dkf2(k-kQ)  =  f  dkN(kQ,a)  -  1. 

—CO 

The  notation  N(kg,o)  stands  for  a  Gaussian  with  mean  at  k^  and  variance 


The  Gaussian  free  particle  wave  packet  is 

-iiSs2  t 

Vx,t)./^(Mo)e-“-e  2" 

The  behavior  of  such  a  wavepacket  is  well  known;  at  t*0  it  i.as  the 


form 


^cj(x*0) 


(8ra2)1/4 


2  2 

-ikQx  e-a  x 


and  the  time-development  of  the  wave  packet  is  well  described  by  the 


From  these  expressions  one  concludes  that  the  centroid  of  a 
gaussian  wavepacket  moves  according  to  Newton's  law  and  is  unaffected 
by  the  momentum  uncertainty,  and  that  the  uncertainty  of  its  position 
is  inversely  proportional  to  the  momentum  uncertainty.  Moreover,  since 
this  observation  holds  no  matter  how  narrow  the  packet,  it  is  not 
implausible  to  pass  to  the  limit  of  an  infinitely  narrow  packet,  that 
is  to  say,  to  a  momentum  eigenstate  |kQ>,  and  to  attribute  to  a  parti¬ 
cle  in  such  a  state  a  quasi-newtonian  motion,  albeit  with  an  infinitt 
position  uncertainty. 

This  argument  can  be  made  more  formally,  and  we  will  now  shew  that 
it  is  possible  to  construct  physically  meaningful  wavepackets  which  are 
equivalent  to  plane  waves,  Bloch  states,  or  Houston  functions,  and 
which  form  sets  of  orthogonal  basis  states.  The  procedure  is  to  gene¬ 
rate  these  states  from  gaussian  wavepackets  whose  variance  approaches 
zero  so  that 

“Vo  lVk-ko>l2  ■  s(k-v- 

To  begin  with,  we  demonstrate  the  orthogonality  of  the  states.  For  any 
wave  packets 


<p  (kQ,x,t)  -  J  dkf(k-kQ)|k> 

—  00 

of  5-function  normalized  states  satisfying  <k' |k>  =  <$(k'-k)  we  have 


/ 


dxij»*  (k^ ,  x ,  t)  V  (kQ  ,x ,  t ) 


■  ^"dk^^dkf’Hk'-kQ^k'  |k>f  (k-kQ) 


•f  dkf *  (k-kp)  f  (k-kQ) 


and  if  we  use  gaussian  wave  packets 


<k-V 


Vk-V  -  (2,)-1/4a-1/2  e  40 


to  form  packets  (kn,x,t),  then 
a  u 


(k-k’)2+(k-k0)2 


dk 


<kQ,x,t|k0>x,t>  =  J"  4/*ip^dx  =  (2tt)  1^2cj 


4  a 


(ko~koy 


(2n)-1/V1  e 


8a 


/ 


2a 


dk  *  e 


8a 


so  that 

lim  J  dx  -  5  .  ,  . 

a+0  J  0  °  kOkO 

This  demonstrates  the  orthogonality  of  the  states.  It  also  pro¬ 
vides  a  basis  for  the  evaluation  of  matrix  elements  in  which  singulari 


ties  can  be  avoided  until  a  final  limiting  procedure. 


Thus  we  confirm  the  validity  of  the  informal  interpretation 
advanced  above  of  the  free  particle  motion  represented  by  a  momentum 
eigenstate.  We  can  provide  a  basis  for  a  similar  treatment  of  the 
motion  of  an  unlocalized  crystal  electron  by  a  similar  development  in 
Bloch  functions  for  the  field-free  crystal,  and  Houston  functions  in 
the  presence  of  a  field. 

—  —  ilc*r  — 

The  results  for  Bloch  functions  >Kk,r)  =  e  u^r)  can  be 
worked  out  using  the  results  of  Appendix  6.3.  We  find 

<V^k0^k^a-*-(/k0^>  =  k0 

(a2+ko)  ■  k0 

a-K) 

<1*'o-*-0^k0^x^a-*-0^k0^>  “  Xnn^k0^  +  h  3k-  ^  ^ 

x 

<  V»0  (k0)  1  x2 1  K+O  (k0}  >_ 1  *  V0  (V  1 X !  (V "  1  2 

3u*(kQ)  3u(kQ)  _ 

3k  3k  dr 

x  x 

and  ft  is  the  volume  of  a  unit  cell. 

The  calculation  for  Houston  functions  wave  packets  differs  from 
that  for  Bloch  functions  only  in  that  the  time  dependence  of  the  Bloch 


where  S(kQ)  =  J 


function 


is  replaced  by 


“■>  [-  feX k*  +  Ea)(k’)dk*] 

*-  x 


so  that  the  differentiation  with  respect  to  kx  yields 


x 


and  we  find,  for  infinitely  narrow  gaussian  wave  packets  starting  from 
rest 

<x>a+0  -  fe  [E  (*?)  -  E<°>]  -  X„n(0) 

<X^> 

cr-H) 

Thus  we  finally  verify  that  in  the  absence  of  collisions  or  tun¬ 
neling  a  Houston  function  does  indeed  describe  Zener  oscillations  with 
a  precisely  defined  phase  and  completely  undeterminate  position. 

We  originally  constructed  the  Houston  representation  as  an 
equal-weight  superposition  of  all  the  Stark  states.  It  is  of  interest 
to  invert  this  procedure  and  determine  the  superposition  of  Houston 


functions  that  is  needed  to  form  a  Stark  state. 


The  Houston  function 


'Kr.t) 


i_  f 

&  *(*£■)*  kJ 


kx(°)  +  ^  E(1)  (k')dk^ 


kx(°) 


can  be  written  in  the  form 


f 

=  »/^/dkx6  ^  kx~kx(°)  ”  ^(k,r)e  ^  *3 


k 


X  E(1)(k') 


(0) 


and  with  the  Fourier  series  representation  for  the  6-function 


( W>  ■f)^E“'1[T(w°>-f)] 


a  wave  packet  of  Houston  functions  with  weight  f(k  (0))  becomes 


i=/dV(k,r) / dkx(0)f(kx(0))E  exp  ^  (kx-kx(0)  -  ^ 


-  —  f°  E(1)at’)dk’  -  —  f  X  E(1)(k’)dk* 

*e  e‘-'kx(0)  *'•'0 


y*dkx^(k,r)  ^  <{> v(k, t )  dkx(0)f (kx(0))e 


i2irv 


ic  x 


k  (0) 


+  -f 

e  *  0 


kx(0)  (1)  - 

x  EU'(k')dk' 
x 


where  we  have  recognized  the  tunctional  form  of  the  Stark  wave  func 
tions  as  in  Equation  (8). 

We  see  that  a  wave  packet  of  Houston  states  ^(kxQ,r,t)  of  the 

form 


/dk  (O)f(k  (O)H(k  n,r,t) 
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is  equivalent  to  a  wave  packet  of  Stark  states  J  dk(J>v(k,tH(k,r)  of 
the  form 


y* 


dk<|)  (k,t)^(k,r) 


—  /*kx(0)  E(1)(k')dk’  -  k  (0) 

ee  Jn  x  *  *  x 


if  g(v)  -  f  dkx(0)f(kx(0))e  ee  *'0 


x  <  x 
e 


the  coefficient  in  the  Fourier  expansion  of 


f(kx(0))exp[^^kx(0)  E(1)(k')dk^J. 


Therefore  to  construct  a  single  Stark  state,  say  with  index  Vq,  energy 
E  ,  we  must  have  g(v)  =  6  which  will  result  if 

V  VV0 


f(kx(0))  =  exp 


Hf 


,(0)  \  _ 

x  Eu;(k’)dk^ 


]exp[.i22™  kx(0)]. 


3.3  STARK  STATES 

We  have  introduced  the  Stark  states  in  Chapter  2,  2.1.3  above,  and 
related  them  to  the  Houston  states  in  that  section,  and  in  Chapter  3, 
3.2.  The  wave  function  of  a  Stark  state  is  localized  by  the  require¬ 
ment  that  the  kinetic  energy  of  a  carrier  in  such  a  state  be 
non-negative.  In  the  present  section  we  shall  study  the  details  of 


this  localization. 


Apart  from  tunneling,  the  Stark  states  are  stationary  states  of 
crystal  electrons  in  an  applied  field.  The  question  whether  a  crystal 
has  exactly  stationary  states  in  a  field  (the  existence  of  "closed 
bands")  has  been  discussed  by  Wannier  and  Fredkin  (9),  but  this  ques¬ 
tion  has  little  bearing  on  the  problem  we  are  addressing.  We  can  view 
tunneling  as  a  perturbation  that  produces  an  energy  uncertainty  in  the 
Stark  levels.  The  magnitude  of  this  uncertainty  affects  the 


observability  of  the  Stark  ladder,  and  it  has  been  argued  (4)  we  believe 
incorrectly,  that  it  is  large  enough  to  destroy  the  ladder.  Indepen¬ 
dently  of  this  argument,  however,  we  shall  show  that  in  the  moderate 
fields  of  interest  to  us,  the  carrier  localization  is  negligibly 
affected  by  the  energy  uncertainty. 

The  wave  function  of  a  Stark  state  can  be  expressed  as  a  super¬ 
position  of  Bloch  functions 
-  £  E  t 

if»v(r,t)  =»  e  v  ^(r) 

r<n 

where  \ jj  (r)  *  /  dk  $  (k)^(k,r) 

v  •'-k/2  x  v 

and  <f>v(k)  =  <~1/2  exP^“(Evkx  ~  f  x  E(1)(P)dk^ 

As  in  Equation  (7),  we  assume  definite  values  of  k  ,  k  .  Such  a  choice 

y  z 

is  in  no  way  restrictive,  since  k  .  *  (0,k  ,k  )  is  conserved  throughout 

^  y  z 

all  computation. 

The  Stark  state  i|>vn(r)  *-s  an  eigenstate  of  the  one-band  Hamilto¬ 
nian  of  the  nth  band 

HJ  (r)  =»  E  ip  (r) . 
n  vn  vn  n 

More  precisely,  this  expression  means  that  the  E^  are  eigenvalues 
of  the  Schroedinger  equation  if  the  interband  terms  are  neglected.  The 
statement  can  be  interpreted  to  mean  that  ,  Ey  are  the  eigenfunctions 
and  eigenvalues  of  a  truncated  Hamiltonian  H  whose  (field-free  Bloch) 

tl 

eigenfunctions  have  X^^O,  m^n.  In  keeping  with  this  interpretation, 
it  can  be  shown  that  the  are  orthogonal  and  form  a  basis  (Appendix 
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6.1).  We  also  show  in  Appendix  6.2  that  the  Stark  levels  E  give  the 


expectation  values  of  the  complete  Hamiltonian 


^’^IHtotal'V1’^  "  Ev 


where  includes  the  untruncated  crystal  Hamiltonian  and  the 


electric  field 


“total  ’  +  V(rt  -  *«• 


The  energy  eigenvalue  of  the  Stark  state  is  given  by 


E  .2^51  + 
V  K 


iy</2  „(!),* 


E'  ; (k)dk 


-k/2 


2™*f  +  <E(1)> 


where  we  have  defined  the  notation 
•k/2 


<f> 


-k/2 


fdk  . 
x 


The  localization  of  the  state  is  described  by  the  expectation 
value  of  the  position  <v|x|v>  and  by  its  root  mean  square  deviation 


Ax  *  [<v|x2jv>  -  |  <v[x|  v>  | 


We  show  this  calculation  in  detail,  since  it  is  useful  in  the 


study  of  Stark  state  wave  packets. 


<v|x  j  v>  *  JJ  dk/dk^^Ck1)  4>v(k)  J  dr 41*  (k'  ,r)x^  (k,r) 


f<nKKHV'>f <n it +  x™(k,)*v® 
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■c/2 
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</2 

<12 
2 


dk  4  *(k)  I  (E  -E(1)(E))  +  (£)|<|,  (k) 

xrv  ee  v  nn  |Tv 


-•c/2 


>]■ 


-  /*  1"-^-  (E  -E(1)(k)  +  X  (k)l  dka 

k  /  ,  ee  v  nn  I  a 

*'-k/2  l  j 


—  +  <X  > 
ic  nn  k 


where  we  have  used  d>  *<P  m  < 

Tv  v 


-1 


Similar ily 


<\j  x  v> 


fK  /dkMk^Mk')^®^*  dr^CiE1  ,r)x%(k,r) 
-</2  J 

r:  + ir+i^r+:n«) 


-</2 


Now 


5T*»(k>  ■[«  (VE<1)(W]  »v(k> 
X 


eE  iS5Tm  *  ® 

x 


o  o 


<v|x  |v>  «  -J_ 


dk 

</2  x 


(E>E-  Ck)--2E- (k)Ev)-  - 
.e  q  x 


+  2iX  j-  (E  -E(1)  (k))  +  i  +  5  (k) 

nn  ee  v  3k  nn 

x  J 

Since  k  is  a  reciprocal  lattice  vector,  the  imaginary  terms  prop¬ 
erly  vanish.  The  expression  can  be  simplified  somewhat  by  using 
E(1)(k)  -  E(k)  -  eeX^,  yielding 


<v  I  *2 1  v>  -  <-fj  (Ev-E(k)  )  2-X^  <k)+=  <k)  >K 

e  e 


To  compute  Ax,  we  write 


<v|x|v>  -  -  —  +  <X  >  - - -  +  —  <E(1)(k)>  +  <X  > 

11  k  im  k  ee  ee  <  nn  < 


-  --  (-E  +<E(k)>  ) 
ee  v  < 

so 

I<v|x|v>|2  -  -j-j  CE2+|<E(k)>J2-2Ev<E(k)><) 
e  e 

and 

<ta>2  '-bh2®’,  -  l<E(k)>J2]+ 

e  E  u  J 

This  result  is  easily  interpreted.  If  we  consider  the  "classical" 
Zener  trajectory,  described  by 

x(t)-x(0)  -  ^  [E(k(t))-E(k(0))] 

where  k(t)  ■  k(0)  +  ^  eet,  we  can  compute  the  mean  square  displacement 


along  this  trajectory  as  a  time  average.  We  find 
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<x(t)> 


1_  r  ee  r* 

ee  IkJ 

L  0 

-\~f 

ee  1iK  J 
L  0 

fe[ 


T  _  _ 

E(k(t) )dt-E(k(0) ) 


] 


dk 


E(kCt))  dkld?  '  E(k(0)) 


] 


K  E(k(t) )dk  -  E(k(0)) 


] 


<E(k)><  -  E(k(0) ) 


] 


where  we  have  used  T  *  77  ,  the  Zener  oscillation  period.  Similarly 


<x2(t)>  -  ~yj  f  fdt  ^ E2 (k) -2E (k) E (0) -E2 (0) 

e  e  J  L 


] 


2  2 
e  e 


*[l 


<  J2.  , 


E  (k)dk 


„  dk  /dt 
0  x 


•/ 


<  E(k)dk 

“2E(0)/  -dWdT+E(0)_ 


2 

e  e 


2 

e  e 


h[  kJ 

hi 


V(k)dk  -2E(0)  7  E(k)dk  +E2(0) 

o  */o 


] 


<E2 (k) >  -  2E(0) <E(k)>  +  E2(0)  1 

K  K  J 


and 


<xZ(t)>  -  I <x(t)> | ^  -  [<E2(k)>  -  2E(0)<E(k)>  +  E2(0) 


-  | <E(k) >  |2  -  E2(0)  +  2E(0)<E(k)>  ] 

^  < 

■  TI  [<e2®>,  '  l<E<k»J2] 

e  e 

which  is  seen  to  be  identical  with  the  squared  uncertainty  of  position 
of  a  Stark  state,  except  for  a  small  term  associated  with  the  polariza¬ 


tion  of  the  band. 
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The  crystal  momentum  expectation  and  uncertainty  in  a  Stark  state 

I  is 

<v|k  lv>  «  0 
1  x 1 


(Ak  )2 
x 


This  indicates  that  the  crystal  momentum  is  entirely  indeterminate 
along  the  electron's  trajectory  in  the  Brillouin  zone. 

We  can  now  see  that  both  one-band  Houston  and  Stark  wave  functions 
describe  the  Zener  oscillations  of  a  crystal  electron,  but  in  rather 
different  ways.  In  Houston  states  the  crystal  momentum,  which  repre¬ 
sents  the  phase  of  the  Zener  oscillation,  is  sharply  defined,  at  the 
cost  of  complete  delocalization  of  the  electron.  In  a  Stark  state,  the 
energy  Is  sharply  defined,  and  the  electron  is  localized  to  the  extent 
determined  by  the  requirement  of  a  positive  kinetic  energy,  but  the 
phase  of  the  oscillation  is  random.  According  to  Equation  (11) ,  the 
electron  oscillation  is  centered  about  a  lattice  site,  with  a  small 
shift  due  to  the  band's  polarization.  The  amplitude  of  the  oscilla¬ 
tion,  for  the  conduction  band  of  GaAs,  is  of  the  order  of  200  lattice 
parameters  in  a  field  of  200  kv/cm. 

3.4  WAVE  PACKETS  OF  STARK  STATES 

The  most  appropriate  representation  of  a  quantum  mechanical  system 
is  normally  determined  by  the  experiment  that  is  to  be  described.  We 
have  already  suggested  earlier  that  an  electron  thermally  excited 
across  the  bandgap  is  unlocalized  and  might  be  represented  by  a  Houston 
function.  The  excitation  will  usually  be  to  the  edge  of  a  band,  k  *  0 
in  a  direct  band  gap  semiconductor.  This  fixes  the  phase,  but  it 


ignores  the  randomness  of  the  excitation.  One  way  to  prepare  a  Stark 
state  might  be  by  Auger  tunneling  from  a  bound  impurity  level,  giving  a 
precisely  defined  total  energy,  again  with  a  random  phase. 

Although  both  of  these  states  represent  electrons  executing  Zener 
oscillations,  evidently  neither  is  suited  for  the  generation  of  coherent 
radiation.*  As  we  pointed  out  in  Chapter  3,  3.1,  we  can  hope  to  find 
a  coherent  superposition  of  states  in  which  both  the  position  and  the 
crystal  momentum  of  the  electrons  are  specified  within  the  limits  per¬ 
mitted  by  the  uncertainty  principle.  Leaving  aside  for  the  time  being 
the  engineering  problem  of  how  such  a  state  is  to  be  prepared  experi¬ 
mentally,  we  now  discuss  a  possible  way  of  constructing  it  from  Stark 
states. 

3.4.1  Minimum  Uncertainty  Product 

For  a  pair  of  operators  A,  B,  with  commutator 

[A,B]  -  iK 

the  uncertainty  relation  is 
AAAB  >  |  K  . 

We  wish  to  construct  a  minimum  uncertainty  wave  function,  for  which 
AAAB  "  |k  . 

*An  extreme  (and  rather  ludicrous)  example  of  this  unsuitability  of 
Houston  states  with  a  range  of  phases  is  offered  by  the  filled  valence 
band  of  an  insulator.  In  an  applied  field,  the%motion  of  every  elec¬ 
tron  in  this  band  is  governed  by  the  equation  hk*qe  ,  and  it  executes 
Zener  oscillations.  Furthermore,  none  of  the  electrons  can  be  scat¬ 
tered,  since  there  are  no  empty  final  states  available.  But,  of 
course,  the  oscillations  are  unobservable,  because  there  is  perfect 
phase  cancellation. 


A  gaussian  wave  packet  of  plane  waves  can  describe  the  motion  of 
an  unaccelerated  particle  in  vacuum.  Such  a  wave  packet  can  be  con¬ 


structed  to  have  a  minimum  uncertainty  product  of  position  and  momentum 
at  a  given  time,  but  because  the  Schroedinger  equation  in  vacuum  is  dis¬ 
persive,  the  wave  packet  spreads  and  the  minimum  uncertainty  product 
grows  with  time. 

In  the  presence  of  a  potential,  there  can  exist  minimum  uncertainty 
wave  packets  that  do  not  spread  with  time.  An  example  is  furnished  by 
the  Glauber  states  in  the  quadratic  potential  of  a  harmonic  oscillator. 
As  we  mentioned  in  Chapter  3,  3.1,  these  eigenstates  of  the  destruction 
operator  do  not  spread,  and  they  have  a  minimum  uncertainty  product  of 
position  and  momentum,  as  well  as  of  occupation  number  and  phase.  When 
the  harmonic  oscillator  is  a  mode  of  the  electromagnetic  field,  a 
Glauber  state  corresponds  to  maximally  coherent  radiation,  as  exempli¬ 
fied  by  laser  light. 

The  minimum  uncertainty  product  of  the  position  x  and  crystal 
momentum  kx  can  be  obtained  from  the  commutator;  for  example,  in  a 
momentum  representation 


3T  +  Xnn<k> 
x 


so 


[x,  kx]<Kk)  -  xkx<j>-kxx$ 


i<J>+ik 


+ 

9k 


X  k  d>-ik 
nn  x  a 


3<f> 


-  k  X  * 
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and  therefore 

AxAk  >  \  . 

x—  2 

3.4.2  Periodicity  of  Expectation  Values 

We  propose  to  construct  a  normalized  wave  packet  of  Stark  states 

—  "  —  ~  W  E  t 

Kr,t)  ■  £  f (vQ-v)^^(r)e 

V=— oo 

with  £  |f(vQ-v)j2  -  1 

—CO 

and  with  the  additional  constraint,  imposed  in  order  to  simplify  the 
calculation  of  certain  averages,  that  the  weighting  function  be  symme¬ 
tric  about  Vq, 

f (vQ-v)-f (vq+v) 

We  can  generalize  the  calculation  of  Chapter  3,  3.3  of  the  expec¬ 
tation  values  and  uncertainties  of  position  and  crystal  momentum. 

Let  A  be  an  operator  having  the  property 

<k'  ,r  |  A  j  k,r>d>  (k)  =  6  (k' -k  )  f"G,(k  )<j>  (k)+C  <fi  (k)l 
»  i  i  >  Tv'  x  x  |_  v  x  v  v  v  J 

with  |k,r>  Bloch  functions  and  the  Stark  state  envelopes  as  defined 

2  2 

above.  The  operators  x,x  ,  k  ,  k  are  all  of  this  type.  We  can  now 

X  X 

establish  that  the  expectation  value  of  such  an  operator  for  a  Stark 
state  wave  packet  <Kr,t)  Is  periodic  with  the  period  of  the  Zener 


oscillation: 
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«Kr,t)  |A|i|i(r,t)>  =  <A> 


/*</2  _  _ _  it  (E  ~E  ) t 

dr  ]T  f*(v  -u)  /  dk'^*(kf)^(kTr)A  e  v 
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and  hence  the  last  term  vanishes  except  when  y=v,  and 

<A>f  -rZf*(v^f(vO“v)e*  V  W  *  <  /*dkxGv(kx)e 


i2it  (v-y)— 


v  u  v 
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With  a“v-y  we  can  rewrite  this 


/2-rree 


<A> 


E  e1"'  ^  tEf*(v„-v+a)f(v0-v)  i /</2 


i2-rock 


f  ®  v vg 

v  a  v 


+  E  ln»0-v)|2cvi 


-K/2 


dk  G  (k  )e 

XV  X 


p 


i 

I 

;< 

► 

! 

1 

■  i 

1 

if 

» 

■  i 


-34- 

which  can  be  recognized  as  a  Fourier  series  representing  a  periodic 
function  of  period  -hx/ee  (the  Zener  oscillation  period),  with  coeffi¬ 
cients 


£f*(v0-^)f(v0-v)  [Fa(Gv)+  Cv6aQ] 


where 

k 

i  r<t2  i27TCt 

F  (F(k  ))  =  -  /  dk  F (k  )e  K 

a  x  <  J  .  x  x 
-</2 

This  means  that  for  such  a  wave  packet  of  Stark  states  not  only  the 
expectation  values  <x>  and  <k^>  recur  periodically,  as  might  indeed 
have  been  anticipated,  but  also  all  their  moments.  The  wave  packet 
reforms  to  its  original  shape  after  a  period  of  oscillation.  It  may 
deform  as  it  oscillates,  but  there  is  no  long-term  spreading.  In  this 
respect  the  Zener  oscillation  is  similar  to  the  motion  of  a  harmonic 
oscillator,  and  unlike  that  of  a  free  particle. 

The  reason  for  the  periodicity  in  time  of  the  packet  wave  function 
is  the  uniform  energy  spacing  of  the  Stark  levels.  Since  the  time  de¬ 
pendence  of  each  of  the  wave  functions  in  the  superposition  is  of  the 
form  exp  [(i/ft)E^t]  with  each  energy  an  integral  multiple  of  h  uzener> 
tha  periodicity  follows  immediately. 

3.4.3  Position  and  Momentum  Uncertainties 

The  time-dependent  crystal  momentum  expectation  value  of  k^  is 


<k  >, 
x  f 


g(o)F  (k  ) 

oc  x 


co 

where  g(a)  =  £  f*(vQ-v+a) f (vQ-v) 

V  OO 


(13) 


and  (see  Appendix  6.4) 


» 
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Similarly  (Appendix  6.4) 


<k2>f 
x  f 


K 

12 


+  Z  c-Da 


ict 
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2  2 
2tt a 


( 2^eE  \ 

'  **  '  C  g(a) 


The  spread  in  of  the  wave  packet,  given  by 


(14) 


(Akx(t))2  =  <k2(t)>f  -  |<kx(t)>|2 

v  v 

varies  in  a  complicated  fashion  in  the  course  of  an  oscillation  period, 
but  because  of  the  periodicity  of  the  wave  packet,  it  is  easy  to  obtain 
a  time  average.  We  have 
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so  that 

-  g  -  5  (f,)2  a5, 

0  a^O  a 

To  compute  the  expectation  value  of  x  we  use  the  result  (see  Equation 


(9)) 


<k,r|x|k,r>4>  (k)  -  [e  -  E(k)]d>  (k) 


so  that 
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g(a)Fa(E(k)) 


The  term  in  E  can  be  simplified  if  we  invoke  the  postulated  symmetry 


f(vQ-v)  =  f(vQ+v): 
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time  average  over  one  period  of  |<x>£  [  is 


hi 


]  dt | <xf  > 
0  v 


TT  Ev  +  E  |  g Cot)  [ 2  I F  (E (k)  )  [  2-2Ev  <E(k)>K 
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Finally  the  expectation  value  of  x  is 
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The  first  term  of  this  rather  complicated  expression  can  be  simplified 

■  Evj'(0)i2  +  i :  i£cc.>]2(E2  +  e2 

y  0  a-1  0  0 
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2 

The  time  average  of  <x  (t)>,  is  (see  Equation  (12)) 

r  v 
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The  squared  uncertainty  of  position  is 
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(Ax(t) ) 2  -  <x2>f  “  l<x>f  i2 

V  V 

with  time  average 

dt(Ax(t))2  =  2  f;  (  ~  j2Qt2  |  f  (ct)  |  2  2  £  |g(a)|2F  (E(k))|2 

*'0  a=l  '  '  e  e  ct=-«° 

-  ~h  [<e®2>k1  -  *  h  <k>  (i6> 

eel-  J 

3.4.4  Uncertainty  Product  and  Band  Structure 

We  have  now  worked  out  the  expectation  values  for  wave  packets  of 
2  2 

Stark  states  of  kx>  k^,  x,  and  x  .  These  quantities  are  periodic  func¬ 
tions  of  time,  with  the  Zener  oscillation  period.  If  the  weighting 

function  f(vrt~v)  goes  to  6  ,  the  expectation  values  reduce  to  those 

0  \>\> 

0 

of  a  Stark  state. 

If  the  uncertainty  product  &x(t)Akx(t)  *  I(t)  takes  the  value 
I(tQ),  at  time  tQ,  it  will  in  general  change  with  time,  and  will  return 
to  I(tQ)  at  t  =  tQ  +  T,  where  T  ■  ich/ee.  Thus  if  a  weighting  function 
is  chosen  to  minimize  I(tg) ,  it  will  not  in  general  produce  a 
minimum  uncertainty  product  at  other  times.  Some  optimization  cri¬ 
terion  should  be  chosen.  The  criterion  will  presumably  depend  on  the 
application  to  be  made,  but  some  general  remarks  can  be  made. 

One  may  wish  to  minimize  the  time  average  of  I: 

\J  I(t)dt  *  <I>T 

or  some  combination  of  and  the  fluctuations  of  I(t)  over  a  cycle 

Yy*T(I2(t)  “  <I>2)dt 

or  the  maximum  value  sup  I(t). 


» 
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b 


Alternatively  one  might  wish  to  minimize  the  spread  of  kx  or  x. 


In  every  instance,  the  quantity  to  be  deduced  is  the  weight  function 
f(v).  The  determination  of  this  weight  function  is  now  a  mathematical 
design  problem  related  to  the  optimization  criterion  that  has  been 
chosen. 

In  connection  with  this  problem  it  should  be  observed  that  the 


expectation  values  we  have  computed  contain  the  Fourier  expansions  of 


—  2  — 

the  band  structure  and  polarization  matrix  elements  E(k),  X^n(k)  and 


E(k).  A  band  structure  of  interest,  that  of  the  conduction  band  of 
GaAs,  can  be  quite  adequately  described  by  five  Fourier  coefficients 
(see  Chapter  4  below).  We  also  know,  from  the  discussion  in  Chapters 
2,  2.3  and  3,  3.3,  that  the  polarization  matrix  elements  are  negligible 
compared  with  the  contribution  of  the  band  structure. 

We  can.  rewrite  the  squared  momentum  and  position  uncertainties, 
time  averaged  over  an  oscillation  periods,  from  Equations  (15)  and  (16): 


U 


(Akx(t))dt  -  Y2  - 


<1.  f  LiiCcOL2 
2tt2  tl  a2 


(17) 


\S 


(Ax(t) )  dt 


22[ 


e  e 


<E(k)2>  -  I <E(k)>  1 2  +  <H  -X2  > 

k  1  k  1  J  nn  nn  < 


a= 


-  E  I  8 (a) I  2 1 F  (E(k) ) 1 2 

a*0 

O--00 


+  E  ct2  j  f  (ot)  [  2 

o-l  K 


(18) 


In  these  formulas,  the  function  g(ct),  it  will  be  recalled,  is  defined 
in  terras  of  the  weighting  function  f(vg-v)  as 
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g(a)  “  23  f*(v0-v+ct)f(v0-v)- 

This  can  be  viewed  as  the  autoconvolution  of  f(Vg-v),  and  there¬ 
fore  g(a)  will  be  somewhat  broader  than  f(Vg-v). 

The  leading  terms  of  Equations  (17)  and  (18)  are  the  squared  un- 
2  2 

certainties  (Ak^)  and  (Ax)  of  a  Stark  state 
2 

<4Vv  -  12 

«*>*  -  Tpr  [<E®S  -  l<E®^!2]+  <h„„  - 

As  we  have  pointed  out  in  Chapter  3,  3.3,  the  position  uncertainty 
of  a  pure  Stark  level  is  effectively  accounted  for  by  the  size  of  the 
Zener  oscillation  orbit.  We  now  see  from  Equation  (18),  that  forming 
a  wave  packet  of  Stark  states  has  two  consequences:  It  broadens  the 
uncertainty  slightly  (last  term  in  Equation  (18)),  because  adjacent 
"Stark  orbits"  are  displaced  from  each  other  by  (roughly)  a  lattice 
parameter;  and  it  can  lead  to  a  much  larger  reduction  due  to  the  corre¬ 
lation  of  the  orbits,  with  each  harmonic  of  the  band  structure  contri¬ 
buting  to  the  reduction. 

For  example,  for  a  sinusoidal  band,  as  would  be  exhibited  by  an 
extreme  tight-binding  model 

E(k  )  ■  -A  cos  k  a 
x  x 

2 

the  band-structure  dependent  dominant  part  of  (Ax)  (neglecting  the 
small  broadening  due  to  the  last  term  in  Equation  (18))  becomes 


I 


(Ax)2  *  j  <E%  -  <E>^  -  L  |g(o)|z|F  (E)|^  I 

e  e  L 

--TT  (I'!  ls<»l2)  • 

e  e  '  ' 

Furthermore,  we  see  from  Equations  (13),  (14)  and  (17)  that  the  band 
structure  does  not  affect  the  spread  of  crystal  momentum. 

3.4.5  An  Example  in  the  Tight-Binding  Approximation 

The  detailed  analysis  and  design  of  a  wavepacket  representing  an 
experimental  situation  depends  on  the  band  structure  of  the  semiconduc¬ 
tor,  and  will  require  numerical  work.  Still  it  should  be  possible  to 
obtain  insights  and  observe  trends  which  might  be  independent  of  the 
particular  material,  using  simple  models  and  analytical  methods.  One 
such  model,  the  sinusoidal  tight-binding  band  structure,  was  introduced 
in  the  last  paragraph.  We  now  study  an  example  of  the  behavior  of 
wave  packets  in  a  model  solid  with  this  band  structure,  a  Real  Equal 
Weight  Packet. 

We  assume  a  superposition  of  Stark  states  centered  on  Vq  with 
weights. 


-(N-1)/2<(v-vq)<(N-1)/2 


f(v-vQ) 


otherwise 


It  is  easily  found  that 


-42- 


g(a)  *  £  f*(vQ-v+a)f (vQ-v)  -  (N-a)/N 
v 

F±1[E(k)]  *  A/2  all  other  F^ [e (k)]  =  0 

FQ[E2(k)]  -  A2/2,  F+2[E2(k)]  »  A2/4,  all  other  F^E^k)]  =  0 

In  addition  the  following  needed  expression  can  be  reduced  to  a  simpler 
form 

C1  -  2  E  +  Ev+1)f*O0-v-l>f(v0-v) 
v 

- 2  E  +  »«<o 

-  AEv0  53  f * (ct  +  l)f(a) 

a 

-2AE  53  (2a  +  l)f*(a  +  l)f(a) 
a 

-  ^vQ  53  f*(a  +  l)f(a)  «  4Ey  g(l) 

a  ® 

where  we  have  used  f*(-a)  ■  f(a),  true  for  a  real  wavepacket  . 

Let  us  use  these  results  in  order  to  compute  the  expectation 
<z*f^  and  uncertainty  Ax  for  the  rectangular  wavepacket.  We  find  for 
our  tight  binding  band 


I 


if  [-\  +AS(1)  cos  (ip)] 

(  “V(1>  -  I  C1  )  cos 


Applying  these  formulas  to  our  real  rectangular  wavepacket  we  find 


ir  K  ♦  A  (¥)  -  (“?)] 


[  (N-D/2 

2(“*>  {fE 

a-1 

'(i-ilWl 

2/l  (N-2)  1  /N-l\2\  2aeet| 

1  \i  n  - 2  ;  ;cos-T-J 


„  (N-D/2  / 

with  the  series  —  £  a2  -  j 

a-1  '  ° 


-(2 ((N-D/2)2  +  3  +  1 


2  2 
IT  .N  A1  N*  „  . 

12  +  4  +  6  12  for  N  large 


Note  that  the  amplitude  of  the  oscillation  increases  rapidly  with  the 
number  N  of  states  in  the  packet.  The  "centroid"  of  a  single  state  is 


1 
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fixed;  as  more  states  are  superposed,  it  oscillates,  reaching  a  maximum 

amplitude  corresponding  to  the  width  of  the  "classical"  Zener  orbit. 

2 

The  position  uncertainty  Ax  contains  a  band-structure  dependent 

term 

A2  [7  1  1  |n-1\2|Jin-2  1  /  N-l  \  2  \  2aeet  ] 
e2£2  2  "  2  [  N  /  j  +  \  2  N  "  2  (  N  j  j  C0S  h  J 

which  narrows  as  the  number  of  states  is  increased,  and  a  second  term 

2 

which  asymptotically  increases  as  N  for  N  very  large. 

This  is  indeed  what  one  might  expect  on  the  basis  of  plausible 
reasoning.  A  single  Stark  state  represents  an  electron  which  has  a 
time-independent  probability  of  being  found  at  each  point  of  its 
orbit — it  is  "spread  out"  over  the  entire  orbit.  As  other  states  are 
superposed  to  form  something  like  a  coherent  state,  they  interfere 
constructively  in  the  vicinity  of  some  point  along  the  orbit,  the 
probability  density  becomes  localized  and  time  dependent,  more  or  less 
describing  a  classical  oscillating  particle. 

This  localization  can  narrow  and  persist  only  if  the  wave  functions 
of  the  Stark  states  forming  the  packet  overlap  in  space;  that  is  to 
say,  if  the  range  of  Stark  energies  is  less  than  the  energy  width 
2A  of  the  band,  or  the  number  of  states  is  less  than  2A/eea.  Stark 
states  which  are  separated  in  energy  by  more  than  2A  do  not  overlap  and 
cannot  interfere;  they  contribute  to  a  probability  density  extending 
beyond  the  "classical"  Zener  orbit,  and  contribute  to  an  increase  in 
Ax,  which  asymptotically  becomes  proportional  to  the  size  of  the  class¬ 
ically  permitted  region.  The  narrowest  localization  occurs  for 

N  -  1.14  (  — )  2/3. 

\eea  / 


3.4.6  Minimum  Uncertainty  Product  Wavepackets 

The  rectangular  wavepacket  used  in  the  preceding  section  was  chosen 
arbitrarily.  We  should  like  now  to  attempt  to  select  a  wavepacket  to 
minimize  the  uncertainty  product  AxAkx»  To  facilitate  this  attempt  we 
will  first  investigate  the  somewhat  simpler  problem  of  minimum  uncer¬ 
tainty  wavepackets  in  the  field-free  crystal,  in  a  basis  of  Bloch  func¬ 
tions.  Since  a  minimum  uncertainty  product  packet  of  plane  wave  states 
is  known  to  be  strongly  time-dependent,  we  will  only  undertake  to  con¬ 
struct  the  packet  at  a  time  t  at  which  the  phase 

exp  [-(i/h  )(t'-t)E(k)]  of  the  Bloch  wave  is  unity  and  the  wave  packet 
has  the  form 

1  f  </2  - 

■  -  /  dk  4>(k  )i|»(k,r) 

K  -</2  X  X 

We  will  try  now  to  derive  the  required  wavepacket  shape  for  which 
the  minimum  uncertainty  possible  allowed  by  the  Heisenberg  relation,  is 
reached  at  a  time  t. 

In  order  to  attain  this  purpose  we  will  use  a  method  similar  to 
that  used  to  construct  the  Glauber  States  and  the  minimum  uncertainty 
wavepacket  of  plane  waves. 

We  showed  that  in  the  Bloch  representation  £r,k]*  which  implies 
that  D(ri)D(ki)  >_  ~  . 

Let  us  call  $_(k  )  the  state  which  satisfies  the  equality  in  this 

u  X 

relation. 

Let  us  look  at  the  origin  of  the  uncertainty  relation 
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\  -  |  l<*l[x,kJU>l 

■  J  |  <4>  |  [Ax,Akx]|  <f>>  | 


where  Ax  *  x-<<|> |x| <j>> 


Ak  ■  k  -<4>  Ik  |  (j>> 
x  x  1  x1 r 

Using  the  triangle  inequality  we  find 

\m  \  I  <<t>  I  C^x»  Akx]|  4>>  |f_|  <(£  I  AxAk^  |  $>  | 


Using  the  Schwartz  inequality  we  can  further  write: 

1/2 


\  1  |  Ax2|  4>>  £«(>  |  Ak^>  j 


|  <  D(x)D(kx) 


The  state  <j>  which  achieves  the  Schwartz  equality  should  satisfy 
G 


Akxl<,,G(kx)>  "  lb^xl't'r.(kv)> 


G'  x' 


(19) 


where  ib  is  an  arbitrary  imaginary  constant.  We  will  have  to  determine 
which  number  b  realizes  the  triangle  equality. 

We  can  rewrite  equation  (19)  as  follows: 


(kx-ibx)$G(kx)  -  (<$G|  kxi  4>G>~ib < <>G  |  x  1 4>G>  ) q>G  (2° 

We  know  that  if  d>„(k  )  is  also  a  solution  of  the  crystal  Uamiltonian 
G  X 

perturbed  by  the  electric  field,  then  <k  >  and  <x>  will  be  periodic 
functions  of  time.  In  the  ideal  case  4>g(kx)  will  satisfy  the  equality 
at  all  times.  In  a  more  restrictive  case  4>G(kx)  will  only  satisfy  the 
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equality  at  a  given  time.  In  any  case  the  function  4>G(kx),  at  all  times 
or  at  a  given  time,  will  be  a  solution  of  equation  (20). 


Let  us  first  solve  the  simpler  equation  (20)  where 


<k^>  =  K(t)  and  <x>  *  R(t) 


Using  x  *  i  +  Xn^(k)  we  can  rewrite  (20)  as  follows 
x 


( \ +  b  sr  * 


ibX  (k)U_(k  )  =  A ( t) 4>  (k  ) 
nn  /  b  x  x 


with  A(t)  *  K(t)-ibR(t) 


The  general  solution  is 


$  (k  )  -  C(t)e 

u  X 


-  k2  +  k 
2b  x  b 


•'A 


x  X(k)dk 


We  will  now  proceed  to  a  study  of  the  proposed  solution  (22)  of  equation 


(20). 


To  normalize  the  wave  packet,  write 


w 


k  Kk 
+  _J 
2b  b 


k 


+  -r2-  +  ik  R  .  f  x  5 
b  x  i  / 

e  Jn 


X(k)dk 


A(t)  K-ibR  _  K 
using  —  *  — g—  =*  b  -  iR- 


Since  X  (k)  is  real 
nn 


k2 


$*(kx)^(kx)  -  |C|2  e  2b  e 


and  the  normalization  Is  determined  by 


i 


k2ct)  /*KCt)+lc/2  (kx-K(t))‘ 

2b/  . .  . 


1  “  I  cl  2  e  2b 


K(t)-</2 


We  find  immediately 


*  K(t) 


d2(V  =  i  b*  b  -  0 


To  determine  other  expectation  values,  we  rewrite 


w 


C(t)e 


k2  +  Dk  +  i. 
2b  x  x  J 


XX(k,)dk* 


with 


D .  m. ,  MSI .  1R 


Let  us  derive  the  average  <x>. 


«t.GkUG>  -  <oG|i  x(k)|^G> 


•K(t)+</2 


K(t)-</2 


+  LX(k)  +  X(k)  dk 


-  i  £  +  i 

D 


(!-«) 


We  see  that  the  proposed  solution  (22)  is  a  solution  of  equation 
(20)  since  <kx>G  ■  K(t)  and  <x>G  =  R(t).  We  have  tested  ^qC^)  against 
equation  (20).  However,  we  need  to  check  also  the  triangle  inequality. 
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For  this  purpose  let  us  derive  the  deviation  D(x)  of  our  coherent  state 
(see  Appendix  6.3). 

,2 


<4>g|x2|4>g>  *  «j>G|-  ~2~  +  2iX(k)  g|-  -  E (k)  |  <J»^> 

3kx 


Let  us  first  deal  with  < - 2~  >G 

3lc  * 
x 


<— 


,.2  G 

3kx  */K-</2 


i 


K+</2 


♦s<v 


J.  «  +  IX(k))*0(kx)] 


dk 


<— 


3k 


2  G 


“-<[-s+i3rx(k)+(‘ 


X(k)  +(  -  ^  +  D  +  iX(k)  i  j  >G 


■)■] 


Let  us  expand  the  square  term 


_/i+D+lx\2,-^|-D2  +X2(k)+  ^kx  -  2iX(k)  (  -  dJ 


then  we  can  write 


Jv-e  2  V  b  ' 


3kx  G 


/•  K-bc / 2  /  k  \ 

A  /2  ^(k)  "  21x00  \  b  +  D  j  " 


i  j—  X(k)dkx  (23) 
x 


Vfe  can  evaluate  the  first  term  of  the  right  hand  side  of  (23) 

2  2 
k  <tc  > 

.1  x  _2  ,  2D  ,  „  1  ^2  ,  2D  „  x  G 

b  "  b2  _  D  +  b  x  G  ’  b  ”  D  +  bK"  fe2 


I.g.R2.21Rf)+aK(K_iR). 


<k2>r 
x  G 
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1 

b 


D(k  y 


+  R(t)‘ 


f 


Then  <x  >, 


+  21«k>  aT+  =<_k»c 

3k  _  x 


£  -  \  1)2  <kx>  +  +  <H(k) 

b 


x2(k)> 


Recall  that  <>  means 
G 


<f> 


WC+k/2 

*&-k/2 


♦S<Wk*,£(Vdk*  • 


Since  and  K  are  time-dependent,  <>  is  also  time  dependent.  Also 


note  that  <>  4  <> 

G  < 


The  deviation  D  (x)  is 


■  ’ 

* 


D2(x)  *  <x2>g  -  <x>2  *  £  -  ^  °2(kx)  +  <H(k)  ~  X20t)>G 

b 

and  the  uncertainty  product  is 

D2(k  )D2(x)  -  i  D2 (k  )  -  —  D4(k  )  +  D2(k  )<H(k)  -  X2(k)>. 

X  D  X  .  <-  XX  Ijf 

D 

The  sum  of  the  first  two  terms,  and  the  last  term  are  each  positive. 
2 

D(k  )  is  a  function  of  b 
x 


D(kx) 


■c/2  -  -  kVb 

k  e  x  dk 
2  •'-k/2  x  x 


/ 


ic /2 

•'-k/2 


c/2  -k  /b 


dk 


o  if  -  d2<v  i»!(y  -  TJ 


tr*0 


b-*-® 


D(kx)2  D(kx)2 

Some  values  of  — r - ~ —  ■  B(b)  are  shown  in  this  table 

l  b  b2 


b 

2 

.02k 

.08k2 

.18k2 

.32k2 

.  5<2 

B(b) 

.2496 

.2474 

.2162 

.1666 

.1250 

* 


a 


* 


2 

Since  the  uncertainty  I  is  greater  than  1/4  we  deduce  that  the  positive 
term 

D(ky)2<S(k)  -  X2(k)>b  >_a2D(kx)2 

is  greater  than  -  B(b). 

2 

Since  the  term  <E-X  >  depends  on  the  form  of  the  Bloch  functions 

b 

2  1 

it  remains  crystal  dependent  and  is  not  likely  to  yield  I  ■  •£■  for 
b>0.  The  Bloch  state  appears  then  as  the  only  minimum  uncertainty 
state  since  l2(b=0)  ■  -  — b— 2^  *  t-  for  b-*-0. 

D 

Mathematically  the  reason  arises  from  the  fact  that  only  the 

Schwartz  equality  is  always  satisfied  whereas  the  triangle  inequality 

is  only  satisfied  exactly  for  b*0.  For  small  b  these  states  are, 

2  1 

however,  very  close  to  achieving  I  *  .  Let  us  call  them  quasicoher- 

ent  and  let  us  denote  them  |y>b  . 

Does  the  family  jy>  constitute  a  family  of  states  which  minimize 

2 

I  ?  In  the  procedure  chosen  we  know  that  they  are  uniquely  determined. 

2 

However,  there  do  exist  other  procedures  to  minimize  I  .  For  example 


we  could  have  tried  to  satisfy  the  triangle  inequality  first.  We  have 
shown  only  that  the  only  exact  minimum  uncertainty  states  of  this  type 
are  the  Bloch  functions  (b*0) .  The  family  |y>,  are  not  necessarily 


those  states  which  minimize  I  .  The  family  |y>^  is  derived  from  a 

2 

special  criterion.  This  criterion  is:  we  will  try  to  minimize  I  by 
satisfying  the  Schwartz  equality.  It  is  difficult  to  analyze  the  use¬ 
fulness  of  this  criterion.  However,  it  is  possible  to  analyze  the 
result  which  is: 

D(k  )2  D(k  )4 

I - g* - -  +  D(k^)  < z  (k)  -  r(fc)>G 


We  have  already  pointed  out  that  the  family  |y>  cottas  close  to  satisfy- 

b 

2  1  k  2  b 

ing  I  «  —  for  b  small  relative  to  ic(b  <  -rrr  ),  since  D(k  )  -  v  and  the 

4  10  x  2 

polarization  terms  are  weighted  by  D(kx>. 

Since  the  ultimate  purpose  of  this  discussion  is  to  explore  methods 
of  localizing  crystal  electrons  in  the  presence  of  a  field  in  regions  of 
reciprocal  and  direct  space  small  compared  with  the  Brillouin  zone  and 
with  the  size  of  a  Zener  orbit  respectively,  the  |y>^  states  evidently 
provide  a  fruitful  basis  for  the  discussion. 

3.4.7  Gaussian  Wavepackets 

We  now  attempt  to  apply  the  calculation  of  the  last  section  to  the 
complete  crystal-plus-field  Hamiltonian.  Since  the  Stark  states  are  a 
complete  basis  for  the  direction  of  a  reciprocal  lattice  vector,  it 
should  be  possible  to  construct  a  wavepacket  with  the  properties  of  a 
Iy*  state  from  Stark  states,  and  we  shall  attempt  to  approach  this 
criterion  as  closely  as  possible. 


We  begin  with 


*G(t’r) 


£  fG00*v<t,r) 


-  i  /*  </2+K _ “i-T 

23  fr(v)  ~p  /  'Hk,r)e 

G  &  A-<n 


,-i  f  *  E(1)(E)dk  1  Ej  ee 

i  ®eyn  x  e 


^x 


K+</2 


-- 

(k,r)e  ee  *^0 


X  P(l) 


Ew(k) 


£  fG("> 


*{---) 

\\e£  h  / 


where  E  -  — +  <E(1)(k)> 

UK 


(k)>  We  can  write  the  waveform  <j>_(k  ,t)  in 
K  OX 


the  Bloch  representation  as: 


-~f 

t)  -  ~  e  ee«'o 


X  E(1)(k)dk  i<E(1)(k) 


i(kxt)  -  ~ 
G  x 


Z^Sc  _  t  \ 
k  \  ee  "h  / 


+  tv  L  1SL  +  k  ) 
)e  <  \  ^  */ 


£  f(v)e  K  ‘ 
v-— G 

.  2tr 

co  iv  —  s 

Let  us  define  h^(s)  »  7!  e  *  f^(v) 

V*— oo 


with  s  ■  k  -  t~*  t 
x  H 


we  have  f_(v) 


if 

*/K- 


...  -iv  2tt 
</2+K  —  s 

e  hr(s)ds 

<12 


As  we  clearly  see  that  the  term  h(s)  will  never  allow  the  destruc 


tion  of  the  phase  modulation  e 


ezJn 


E(k’)dk' 


if  f  (v)  is  time 


independent  and  hence  <J>  (k  ,t)  will  never  have  the  form  of  the  solution 

x 

(22)  of  our  quasicoherent  state  Iy*,  at  all  times. 

b  2  /  s2 

-a  (v-vQ) 

However,  using  a  gaussian  distribution  e  we  can 

generate  a  periodic  gaussian  function  h-,(s)  centered  at  zero  with  a 

VJ 

small  deviation  such  that  for  s  |k  -  y  ,  y  +  kJ 


/■*  i  v 

J  6 


x  fg(v)dv 


iv  —  s  -o^(v-vn)^ 
e  '  e  0  dv 


•  i  f-  S(V+V  )  _  2  y 

e  e  dV 


.  2rr  -  ~ 

i  —  svn  -2 

k  0  2o 

e  e 


for  K  -  |i  s  <  |  +  K 


we  see  that  h(s)  is  centered  at  zero.  We  also  need  a  to  be  small. 
Let  us  substitute  for  h(s)  in  equation  (24) 


-  ^  fX  E(1)aOdk  i<E(1)(k)> 
*G(kxt)<X  *  J°  e 


te-t) 


K 

,  2tt  /,  est  \ - — s- 

—  v0  lkx  '  — )  2a2 


2  2  2  ectk 

e  E  t  _ Jt 

2  2  2 
h  2a  Ha 

e  e 


except  for  the  replacement  of  Xnn(k)  by  E^(k)  equation  (25)  can  be 
developed  in  the  shape  of  equation  (22) 


with  CCt)®  e 


b  =  cr 


n  *  u  u  ii  ^  u 

e  e 


/-is  _  k  2ttV  eetk 

k  -  i<E(1)(k)>  ~7+i— - -k  + - 5s 

bx  k  ec  kXa.2 


ho 


with  K(t)  *  +  ~ 


1  (1)  — 

R(t)  -  -  —  <Eu'(k)> - - 

ee  k  tc 


we  can  write  ^,(k  )  as 


<j>G(kx)  *  C(t)e 


i  k2  +  m.  k 

2b  x  b  x 


j-jy  E<i><k')dki 


What  are  the  properties  of  this  state?  The  expectation  values  have  been 
calculated  before.  It  is  possible  to  replace  X^Ck)  by  E^(k)  in  the 
average  before  any  summation  is  done  over  (k) .  We  directly  deduce, 
as  before 


"Vg 


K(t) 


D(k  )  2 


2 

CT 


2 

if  c  small 


<x  >  -  ± - f— +  R(t)*  +  <E(k)  -  X  (k)>G 


2R(t)  <g/£\>  +  <  £  (k),  > 

+  e£  +  22  >Q 

e  e 


the  deviation  is  finally 

(  \2  I  °(Vb  2 

D(x)  "  b  "  b2  ~  “  X  Gb(t) 


k2  [ 


2  — 


<E‘(k)  „  -  (<Effi> 


Vc) 


Gb(t))2] 


We  directly  deduce  that  the  Houston  state  is  a  minimum  uncertainty  state 


since 


D(kx)2D(x)2 

b-K) 


f  +  D(kx)(...) 
b-K) 


1 

4 


As  we  have  done  earlier  with  the  rectangular  wave  packet  of  Stark 
states,  we  shall  now  study  the  behavior  of  a  wave  packet  of  Stark  states 
with  gaussian  weighting,  in  a  band  with  the  tight-binding  cosine  form. 
Before  we  proceed  to  details,  we  should  point  out  that  our  calculation 
so  far  has  included  an  arbitrary  restriction  to  wave  packets  with  real 
weighting  coefficients.  As  it  turns  out,  this  restriction  is  not 
altogether  trivial.  For  example,  we  have  seen  in  Equation  (8)  that  in 
order  to  form  a  Houston  state  as  a  superposition  of  Stark  states,  it  is 
necessary  to  use  essentially  complex  expansion  coefficients,  that  is  to 
say,  coefficients  whose  relative  phase  does  not  vanish  at  any  time 
during  the  Zener  cycle. 
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A  consequence  of  the  restriction  in  the  present  context  is  that 
the  minimum  "size  of  the  electron"  Ax  to  which  a  real-coefficient  wave 
packet  of  Stark  states  can  be  localized  is  field-dependent,  and  as  we 
shall  see  presently,  may  be  relatively  large;  while  a  complex  coeffi¬ 
cient  wavepacket  is  only  band-structure  limited,  and  can  be  of  the  size 
of  a  unit  cell  (e.g.,  a  Wannier  function). 

As  a  practical  matter,  however,  the  restriction  is  not  very  serious. 
A  strongly  "compressed"  electron,  as  one  might  expect,  fluctuates  dra¬ 
matically  in  size  over  a  Zener  period,  and  is  far  from  our  notion  of  a 
coherent  state.  At  the  same  time  we  shall  see  that  states  approaching 
this  notion,  whether  formed  with  real  or  complex  weights,  have  very 
similar  properties. 

We  have  determined  that  for  a  real  gaussian  weighting  function 
with  deviation  a  in.  k-space,  the  uncertainty  Ax(t)  oscillates. 

We  shall  now  determine  the  deviation  cr  which  minimizes  both  the  mini¬ 
mum  and  maximum  of  Ax(t) ;  we  already  know  that  for  narrow  wave  packets, 
a  10  \  such  gaussian  wave  packets  come  close  to  being  coherent 

states  with  uncertainty  product  AxAkx  3  0.5. 

From 

Ax2  =  +  -i-j  j  <E2  (ic)  >G(c)  |  <E(k)  >G(t)  |  2  J  for  small  b 

with 

E(k)  *  -A  cos  k  a 
x 

-I  |k  .  eet\  2 

G(t)  ■  - - —  e  k  '  x  ^  for  small  b 


its  variance  is  Akx*  a  =>  /-j  ,  let  K  *  — 


<E(k)> 


G(t) 


1  2 

"  (VK) 

—  -  f  "  e  2a  cos  k  a  dk 


let  y  -  kx-K 


1  2 


<E(k)> 


ZA_I  /* 


co  2ct 

e  cos(y+K)a  dy 


cos  (y+K)a  =  cos  ya  .  cos  Ka  -  sin  ya  •  sin  Ka 


<E(k)>G^t^  »  (-A)  cos  Ka<cos  ya>N^Q  a)-sin  Ka<sin  ya> 


(-A)  cos  Ka  e 


CT2  2 
2"a 


where  we  use  the  expression  of  <cos  ya>N(g  and  <sin  ya> 
tabulated  in  the  appendix  6.5. 

Similarly 


<E2(k)>G(t)  ■  *2<  7  008  2kxa  +  7 


2  2 
A  .  A 


—  +  2“  cos  2Ka  <cos  2ya>N(0,o) 


2  2  2  2 
A  .  A  __  -2c  a^ 

Y~  +  —  cos  2Ka  e 


and  finally 
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.2  1.1  I  A2  ,  A2  / 2eet  \ 

**  +  r  +  rcoshra 

eel 


2 

-2a  a 


.2  2  /eet  \  -a  a 

-  A  cos  I  — : —  ale 


(¥•)■ 


.  2  1  .  A' 

Ax  -  orr  + 


2b  2e2£2 


1  -  e 


ba 

2 


) 


/  2eeta  \  /  -ba2  "  2 

+  C08ph~]  ^ 


The  position  of  this  electron  is  given  by 


fe[< 


<E(k»G(t)  -  Ev] 


b  2 

A  a  eet 

A  e  cos 


2eeta 


The  minimum  value  of  Ax  is  attained  for  cos  — t- —  =  1  which  corresponds 


to  the  edge  of  the  B.Z. 

2  1_  A2 

Axmin  2b  9  2  2 

2e  e 


I- 


ba 

,  “  2  .  -ba" 

2  e  +  e 


For  a  given  electric  field  e  the  minimum  size  of  the  electron:  Ax  is 
given  by 


able  -e  I  -  ^  with  n  =  — a 

o  10^ 

For  a  typical  crystal  2A  -  2eV, a  -  5A  and  we  have  n  -  0.4  — —  with  e 
expressed  in  kV/cm.  We  can  now  plot  the  results  formed  as  a  function 
of  the  electric  field  in  Figure  1.  We  see  that  Ax  is  large  at  t  =  T/< 
and  the  uncertainty  product  is  large  at  that  time. 

Let  us  also  derive  the  state  where  the  maximum  value  of  Ax  is 
minimized 

Ax  max  *  -sr-  + 


dAx2 

db 


0 


2  2 


n  a 


8 


2 

a 


0 


-ba  .24  4 

e  b  a  -  — j 

n 


-ba' 12  .  2 
e  ba 


Figure  2  shows  a  state  with  less  fluctuations  in  spread  Ax;  its  uncer 
tainty  product  remains  close  to  0.5  at  all  times. 
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Figure  2.  Localization  of  a  quasicoherent  state  of  the  tight-binding 
band  structure:  Ax  minimized. 
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We  have  earlier  introduced  the  Houston  functions  (Equation  (8)  in 
Chapter  2,  2.2.3  and  Chapter  3,  3.2)  as  superpositions  of  Stark  states 
with  complex  weighting 


dkx(0)f(kx(0))e 


i/ 


k  (0)  m  - 
x  Eu;(k')dk'  - 


kx(°) 


which  is  the  Fourier  coefficient  of  f (k  (0))  weighted  by 


i  (°)  (1)  - 

-/  x  E(1)(k')dk^ 


only  if  the  waveform  f(kx(0))  contains  a  phase  modulation 


waveforms  f(v)  generate  the  type  of  oscillations  we  have  discussed  up 
to  now. 


.  The  weighting  function  f(v)  is  real  if  and 


k  (0)  m  _ 
x  E^Ck'Mk' 


which  cancels  that  weighting  f(kx(0)).  Real 


We  will  now  consider  a  waveform  f(kx(0))  without  this  phase  modula¬ 
tion.  Such  wavepackets  generate  a  mode  in  which  the  size  of  the  elec¬ 
tron  Ax  is  not  bounded  by  a  minimum  value  set  by  the  electric  field. 

Let  us  recall  some  previous  results  about  wavepackets  of  Houston  states. 
We  will  set  X(k)  =  0  in  this  calculation.  We  have 


<x>  -  i<f*(k  (0))f ’ (k  (0))> 

X  X  K 

+  ir<{E  r+  k(0))“  E(?(0))Jlf{kx(0))l2>< 

2 

<x  >  -  -<f*(k  (0))f*' (k  (0))> 
x  x  < 

+  < [e* ( k(0)  +  -  E1  (k(0) ) |  |f(k  (0))|2> 


+  ~Yz  <|E(k(°)  +  -  E(k(0))j2|f(kx(0))|2>K 

+  ||  < | E CO)  +  -  E(kC0))jf*(k[(0))f,(kxC0))> 


Q 


+  <■  (k(0)  +  |f(kx(0)) 


1  /*</2 

with  <(  )>  »  -  /  dk  (0)  (  ) 

K  J-</2  x 

Let  us  simplify  these  expressions  by  writing 


E(k(0)  +  — )  »  E(t) 


E(k(0) )  =  E(0) 


f(k  (0))  =  f(k)  =  f  =  |f|  eJ0 

We  first  wish  to  prove  that  it  is  possible  to  localize  an  electron 
in  this  mode.  For  this  purpose  let  us  simply  build  a  Wannier  state 
positioned  at  l  along  the  electric  field.  The  corresponding  wavepacket 


-ik  (0 )l 


f(kx(0))  *  e 


Using  a  cosine  bandstructure 


E(k)  *  -A  cos  ak 


we  quickly  find 


1 


Ax  »  <E(k(t))><  +  ~2  2  ^  1  -  cos  — —  j 

At  times  t  *  m  — 2tt  our  electron  is  localized  to  the  extent  allowed 
aee 

by  E(k)  which  is  of  the  order  of  one  lattice  parameter.  Having  done 
this  the  electron  size  will  oscillate  drastically  to  the  size  of  the 
Zener  oscillation.  The  point  here  is  that  the  localization  at  a  time 
t  is  not  electric  field  limited  but  is  bandstructure  E(k)  limited  in 
this  mode.  However,  its  spread  is  electric  field  dependent.  Let  us 
generalize  these  results.  We  have: 

<x2>  -  -<f*f">K  +  -Jy  <{E(t)  -  E(0)  }2  | f  |2> 

e  e  K 

+  ^  <{E'(t)  -  E'  (0) }  |f  |2>|C 


+  ~  <{E(t)  -  E(0) } 

We  provisionally  drop  the  E(k)term 

using  <{  E*(t)  -  E'(0)  }  |f|2>  =  {  E(t)  -  E(0)  }  |f|2]*^2 

_i_<{E(t)  _  ECO)  >  2|f  I  |f  I  »>K 

i 

and  f*f'  -  f*(|f|,e+j<^  +  •  |f \e+i*) 

-  2 1  f  |  (  f  |  T  +  j/|f|2 
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we  find  <Ax2>  ■  -<f*f'’>K 

+  ~J~2  <tE(t)  -  E(0)  )2  |  f  1 2>|C  -  ^  <(E(t)  -  E(0)  }«f> '  |  f  |  2>k 

Since  there  is  no  phase  modulation  in  f,  the  phase  <p  is  only  related  to 
the  initial  position  in  the  crystal:  <f>  »  -£kx(0).  In  the  desired 
term  Ax  all  terms  involving  t  should  cancel  and  we  find  (letting  t  *  0 
for  simplicity) : 

Ax2  -  +  ~YJ  <CE(t)  -  E(0)  )2  |  f  |  2>k 

e  e 

5~2  ( <(E^t)  -  E(0) } | f 1 2><j 2  +  <HCk) |f |2>< 

<x>  -  —  <(E(t)  -  ECO) } | f | 2> 
ee  ic 

2 

For  t  *  0  the  bandstructure  terms  cancel  and  Ax  reduces  to  -<f*fT,>K 
which  is  the  variance  of  the  wavepacket  related  neither  to  the  band- 
structure  nor  to  the  field.  These  last  equations  can  be  compared  with 
those  of  the  superposition  with  real  weights.  In  that  mode  no  E(0) 
terms  appears  as  a  consequence  of  the  phase  modulation 

i  /*k 

■  x  E(k)dkx 

e  .  If  an  electron  Is  not  too  strongly  localized  at  a 

time  t  its  size  will  not  oscillate  so  drastically.  An  interesting 
application  will  be  to  determine  what  is  the  minimum  size  of  an  electron 
traveling  through  the  Brillouin  zone.  In  ballistic  transport  the  de¬ 
vices  considered  are  usually  at  the  order  of  a  micron  or  smaller.  It 
can  be  imagined  that  for  such  device  dimensions  the  size  of  the  electron 
may  not  be  negligible  on  the  scale  of  the  device. 
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The  wavepacket  to  be  used  is  obviously  the  quasicoherent  one 


f(kx(0» 


-(k  (0)-kn)2/4o2  -  i£k  (0) 

X  w  X 


centered  at  kQ  with  Z  *  0  (without  any  restrictions) 


f'(k  (0))  -  f 
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/k  (°)-k  x 


f'»(kx(0))  -  f’ 
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f(kx(0)) 


and  we  find 
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We  look  for  a  wavepacket  of  variance  a  which  will  minimize  Ax 
It  will  turn  out  again  that  such  wavepackets  are  sufficiently  narrow 
for  the  Gaussian  wavepacket  used  to  be  nearly  a  coherent  state 
AxAkx=  0.5.  This  will  assure  our  result  to  be  near  optimum.  In  this 
calculation  we  set  k^  *  0  and  we  use  the  cosine  bandstructure 


E(k)  ■  -A  cos  akx(0)  -  E(Q) 


E(t)  ■  -A  cos  a|kx(0)  +  ^  tj 


For  narrow  wavepackets  the  two  first  terms  of  Ax  reduce  to 

.  2 

Let  us  compute  the  second  term  of  bxT 


<  { E(t)  -  E(0) }  > 


N(0, a) 


<A2{cos  a  (k+K)  -  cos  ak}2>KT/,n  .  with  K  -  ~ 

N(0,  o;  tl 


2  1  1 

*  A  <1  +  —  cos  2a(k+K)  4-  -j  cos  2ak  -  cos  a(2k+K)>-  cos  aK> 


'[1  +  i 


cos  UK  c-*h°2  +  I  -  cos  aK  a'*2  2°2 


-  cos  aK 


where  we  have  used  the  formula  of  appendix  6.5.  The  third  term  of  Ax 
is  the  square  of 


<E(C)  -  E(0)>h(0c)  -  -Aacos  a(fcfK)  -  cos  ak>N(0_o) 


■  -A  e 


(cos  aK  -  1) 


and  we  can  finally  write 


.2  1  .  A 

**  La1  2  2 

4o  e  e 


r  _  2  2  - 

1  -  cos  aK  -  e  a  (1  -  cos  aK) 


2  2  , 

+  e  cos  aK  (cos  aK-l)J  with  K  = 


At  time  t  ■  0 : 


Ax2  "  l/4o2  ■  Ax2  min. 


At  the  Brillouin  zone  edge  t  »  : 
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Ax' 


a2  (  2  2 
_L_  +  _A_  )  ,  _  -  -a  <J 

.  2  +  2  2 
4a  e  e  ( 


,22) 
+  2e-2“  ® 


i_  +  _2ii  I  x  _  -a2b/2  I  2 

2*  .2,2  | 


where  we  Introduce  b  *  2a  or  a 


vr- 


Ax  is  an  interesting  oscillat- 


1  It*  X  7T 

ing  function  of  time.  Its  maximum  occurs  at  t  ■  r  -  ■  r  or  K  *  — 

2  ee  a  4  2a 

which  is  1/4  of  the  Zener  period  T.  This  maximum  is 


2  1  A2  )  .  -a2b/2 

Ax  -  -sr  +  -r-5-  <  1  -  e 


2b 


2  2 


e  e  |  J 

To  find  the  wavepacket  which  minimizes  this  quantity  we  find  the  solu¬ 
tion  of 

¥-0 


1  22  2  2U/0 

or  -  -iy  +  V-  '  I"  • '  0 

2b2  4  2 


-a  b/2  .24  4 

i  b  a  -  -j 
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n 


_2 _ 

-a2b/4 

e 


The  results  found  using  these  formulas  are  the  same  as  for  the 
Zener  mode  when  Ax(T/4)Max  is  minimized,  although  the  expressions  are 
not  exactly  the  same  (the  differences  are  extremely  small).  We  compare 
these  results  in  the  table. 
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Real  Coefficients 


Complex  Coefficients 


Ax  Max  minimized 


Ax  Max  minimized 


by  n 


.  2  -a  b/2 
ba  e 


.  2  -a 
ba  e 


Ax  Max  occurs  at  7- 
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Ax  Max  occurs  at  y 
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Ax  Min  occurs  at  ^ 


—  +  ^fl  - 
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A  2  i_ 

Axmin  "  2b 


Remark 


h  2tt 


*  -  *  Zener  period 


for  n  ■  - .  whic 

.  2  -a2b/4 
ba  e 

minimizes  Ax  Max 


Ax^(T/2) 


1  +  a2/, 

ib +  r|1"€ 


2.  2 

a  b  \ 


therefore 
Ax  Max  -  Ax  Max 
real  complex 
Ax  Min  *  Ax  Min 


This  is  only  true  for  the  Ax  max  minimized  wavepackets! 


The  graph  for  the  dependence  of  Ax  max  minimized  in  the  complex 
coefficient  mode  is  then  the  same  as  graph  2  of  the  real  case.  The 


uncertainty  product  is  nearly  0.5  at  all  times. 

In  order  to  show  that  the  Ax  max  minimized  state  is  a  special  case 
we  have  also  studied  a  different  case  where  the  electron  is  localized 
initially  at  t*0  over  10  lattice  parameters  (Ax  -  10a) .  We  have 
plotted  on  Graph  3  Ax  for  different  times: 


Ax(0)  - 

Ax(10"12 

Ax(T/4) 

Ax(T/2) 


10a  - 

sec) 


Ax 


min 

-12 


10 


sec  -  scattering  time 


■  Ax 

max 

■  Ax  on  the  Brillouin  zone  edges 


On  graph  4  we  show  the  time  dependence  of  Ax  for  Ax(0)  ■  10a  for 
e  *  0.1,  £  ■  1,  e  *  10,  e  =  100  kV/cm. 
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4.  FOURIER  SPECTRUM  OF  THE  GaAs  CONDUCTION  BAND 

4.1  MOTIVATION 

We  have  seen  that  the  coefficients  of  the  Fourier  series  expansion 
(FSC)  of  the  band  structure  in  the  kx  direction  appear  as  parameters  in 
the  general  formula  for  wave  packets  of  Stark  states.  Since  conduction 
electrons  in  GaAs  are  possible  candidates  for  the  observation  of  Zener 
oscillations,  the  FSC  of  the  GaAs  conduction  band  will  play  a  role  in 
the  design  of  wave  packets. 

More  importantly,  these  FSC  are  also  used  for  calculating  the 
radiative  transition  probabilities  between  Scark  levels.  We  shall 
demonstrate  this  by  working  out  the  probability  p^  of  a  transition 
from  an  initial  state  |v>  to  a  final  state  |y>. 

We  assume  the  Stark  states  to  be  exact  eigenstates  of  the  crystal 
plus  electric  field  Hamiltonian  Hq 

H0IV>  =  EJV> 

the  total  Hamiltonian  is  then 


H  „  . 
total 


H, 


+  p»A 


where  A  is  the  vector  potential  of  the  electromagnetic  field  and  p  the 
electron  momentum.  If  we  assume  the  interaction  potv'tial  to  be  suffi¬ 
ciently  small  (compared  to  the  inverse  of  the  observation  time)  for 
first  order  perturbation  theory  to  be  valid,  then  in  the  dipole  approxi¬ 
mation  the  transition  probability  p  is  proportional  to  the  square  of 
the  matrix  element  <y|px|v>.  Let  us  evaluate  this  matrix  element: 


<U|px|u>  -  <ul*Htotal  -  Htotal*|v>  -  (Eu  -  Eu)<U|X|U> 


<Ev  -  V  7  f ,  k  <E(1)(fc)  -  Ev>  +  X„n‘M  dkx 

J  -k/2  L 


(E  -  E  ) 

V  \x 


r.^s  +i  r 

|_  ee  vp  kJ 


f  </2  E (k)  gee  ^Ev  Ep^kx  dfc 

'  ee  x 

-k/2 


cv 

i  1  r<n  i27i(v-y)  -p  _ 

S  -  E  )  —  -  /  e  E (k)  dk 

v  y  ee  k  J  x 


We  conclude  that  this  transition  probability  is  proportional  to 

the  square  of  the  amplitude  of  the  FSC:  F  of  the  bandstructure  E(k). 

v-y 

The  transition  probability  is  the  basis  of  the  phenomenological  inter¬ 
pretation  of  radiation  where  the  ideally  unscattered  Zener  electron 
initially  prepared  in  state  | v>  moves  to  Stark  states  of  lower  energy 
in  the  field  direction  and  releases  its  potential  energy  as  radiation. 

The  process  that  we  have  just  described  is  the  usual  picture  of 
the  incoherent  emission  of  radiation.  For  the  device  application  of 
Zener  oscillations  as  a  coherent  radiation  source  we  are  more  interested 
in  a  classical  description,  in  which  the  conduction  electrons  oscillate 
in  phase.  In  this  classical  picture  the  radiation  originates  from  an 
oscillating  dipole  which  is  formed  from  a  superposition  of  eigenstates. 
We  have  studied  such  superpositions  or  wavepackets  at  length.  The 
additional  feature  is  that  as  it  radiates  the  wavepacket  will  move  and 

K 

spread.  We  know  that  for  a  wavepacket  not  too  wide  in  Akx  <_  the 


expected  position  of  the  electron  is  given  by 
x(t)  -  x(0)  +  —  [E(k(0)  +  t)  -  E(k(0))J  . 
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We  also  found  that  this  electron  was  localized  in  a  size  2Ax  which  is 
also  oscillating.  Most  of  our  studies  emphasized  the  various  kinds  of 
behavior  of  Ax  for  different  wavepackets.  For  the  optimal  wavepacket 
where  AXji^  is  minimized  it  was  seen  that  Ax  does  not  fluctuate  much 
the  product  AxAk^ being  nearly  optimal  at  all  times.  For  such  a  wave- 
packet  the  oscillating  electron  is  a  simple  classical  dipole.  Ax  being 
small  compared  to  the  amplitude  of  oscillation.  The  frequency  of  os¬ 
cillation  of  this  dipole  is  then  given  by  the  Fourier  analysis  of  x(t) 
which  is  simply  the  Fourier  analysis  of  the  bandstructure.  The  Fourier 
analysis  then  yields  the  radiation  spectrum  of  an  unscattered  Zener 
electron.  Our  studies  of  wavepackets  have  enabled  us  to  establish  the 
conditions  under  which  the  unscattered  Zener  electron  exhibits  classical 
behavior . 

4.1.1  Fourier  Series  Computation 

The  bandstructure  of  the  GaAs  conduction  band  was  kindly  furnished 
by  Prof.  Karl  Hess  of  the  University  oi  Illinois.  It  is  given  at  156 
sampled  points  in  (1/48 )th  of  the  GaAs  Brillouin  zone  (see  Figure  5  and 
program  DATABAND) .  For  convenience  we  use  the  bravais  cell  of  the 
reciprocal  space  which  has  simple  boundaries.  Using  the  48  fold  symme¬ 
try  of  the  Brillouin  zone  we  fill  ^  of  this  Bravais  cell  (see  Figure 
6  and  program  FILLBZ2)  and  store  in  the  file  CUBEBAND.  The  Fourier 
computation  is  carried  out  by  the  program  called  ZENER2  (see  Figure  7). 
This  calculation  is  carried  out  for  each  desired  direction  of  periodi¬ 
city  <  in  reciprocal  space  starting  from  any  point  k0  of  the  Brillouin 
Zone.  An  accuracy  of  one  percent  was  required.  This  accuracy  is 
determined  by  the  number  M  of  points  used  for  the  interpolation  of  the 
bandstructure  along  the  path  [<0,  <0  +  < ].  This  interpolation  uses  a 


TYPE:  RUN  ZENER  2 


ZENER2  * 


LOAD  BANDSTRUCTURE  DATA 
FROM  CUBEBAND  FILE 


DIRZENER2 


DIRECTIVE  PROGRAM 
K-  ?  KO=?  M  =  P 


TYPE: RUN 
ZENER2F 


PARA3RD2 


LOAD 

POURFILES 

OR 

BANDFILES 


INTERPOLATION 
2ND  OR  3RD 
ORDER 


DISPATCH2 


ZENER2F 


DISPATCH  PROGRAM 
(MENU) 


FOURIER  SERIES 
FROM  K=  ? 

TO  K=? 


F0URSERIES2 


PLOT  OR  PRINT 
OR  SAVE  INTO  FILES, 
INTERPOLATION  OR 
FOURIER  RESULTS 


BANDPRINT2 

BANPPL0T2 

F0URPRINT2 

F0URPL0T2 

BANBFILE2 

F0URFIIE2 


*  Actual  name  of  the  programs 


Figure  7.  Block  diagram  cf  the  computer  programs  used  to  Fourier 
analyze  the  band  structure  in  an  arbitrary  direction. 


polynomial  expansion  of  3rd  order  using  19  coefficients  fitted  to 
the  local  sampled  data.  The  formula  for  the  point  located  at  (x,y,z) 
of  the  closest  sampled  point  (W1,W2,W3)  is: 

E(x+Wl,  y+W2,  z+W3)  =  E(W1,  W2,  W3) 

+A(0)x  +  A(l)y  +A(2)z 
+B(0)x2  +  B(l)y2  +  B(2)z2 
+G(0)x3  +  G(l)y3  +  G(2)z3 
+F(0)yz  +  F(l)xz  +  F(2)xy 
+H(l)x2y  +  H(2)xy2 
+H(3)yz2  +  H(4)y2z 
+H(5)x2z  +  H(6)xz2 
+H(7)xyz 

The  Fourier  computation  is  made  using  the  trapeze  method.  The  results 
are  saved  in  files  labelled  k,  kO"  and  M.  k  and  <0  are  given  in  Miller 
notation  and  correspond  to  a  body  centered  cubic  reciprocal  lattice. 

In  Figure  7  the  block  diagram  for  these  programs  is  shown.  Those  pro¬ 
grams  are  given  in  appendix  6.6). 

4.1.2  Results 

The  Fourier  calculations  are  given  for  <  ■  (100),  (110),  (111)  on 
pages  82-84  .  Five  harmonics  (10  for  (111))  are  sufficient  for  an 
accurate  description  of  the  GaAs  bandstructure.  In  practice  it  might 


be  interesting  to  use  for  example  the  fourth  harmonic  in  the  (111) 
direction  in  order  to  maximize  the  number  of  oscillations  of  an  electron 
before  it  is  scattered. 

The  same  calculations  were  done  with  at  a  slight  angle  with 
respect  to  a  principal  direction  in  order  to  evaluate  the  effect  of 
misalignment  of  the  field.  As  expected  the  radiation  spectrum  is 
largely  conserved  except  for  a  broadening  and  a  curious  splitting.  It 
is  possible  to  account  for  this  behavior  analytically,  but  a  plausi¬ 
bility  argument  is  a  good  deal  more  instructive. 

If  the  field  is  along  a  principal  lattice  direction,  the  trajectory 
in  the  extended  Brillouin  zone  that  describes  the  motion  of  the  elec¬ 
tron,  beginning  for  example  at  T,  traverses  an  identical  path  in  each 
repetition  of  the  central  zone,  that  is  to  say,  it  is  periodic  with 
period  k,  and  the  Fourier  coefficients  are  determined  by  the  band  struc¬ 
ture  along  this  segment  of  "length"  k. 

If  the  field  deviates  from  a  principal  lattice  direction,  the 
trajectory,  starting  at  r,  may  traverse  many  repetitions  of  the  central 
zone  before  it  again  reaches  the  point  r.  Thus  the  fundamental  period 
of  the  periodic  motion  may  be  long,  to  be  exact,  the  length  of  the 
reciprocal  lattice  vector  of  the  trajectory.  Furthermore,  the  trajec¬ 
tory  will  sample  different  regions  of  the  central  zone  as  it  crosses 
successive  repetitions.  However,  if  the  deviation  from  a  principal 
direction  is  small,  then  because  of  the  continuity  of  the  band  struc¬ 
ture,  successive  segments  of  the  trajectory  will  be  quite  similar, 
changing  only  slowly  in  the  course  of  the  transit  from  r  to  f*  Thus  the 
long-period  periodic  motion  can  be  viewed  equivalently  as  a  short  (ic) 
period  modulated  motion,  with  a  "low-frequency"  modulation 
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corresponding  to  the  gradual  change  of  the  band  structure  segment  that 
is  being  traversed.  From  this  viewpoint,  the  broadening  of  the  spectrum 
arises  from  a  low  frequency  frequency  modulation  of  a  high-frequency 
carrier. 

Furthermore,  although  the  change  is  gradual,  corresponding  to  a 
low  modulation  frequency,  it  can  be  quite  substantial  for  a  band  struc¬ 
ture  as  complicated  as  that  of  the  GaAs  conduction  band:  adjacent 
segments  of  the  trajectory  are  rather  alike,  but  segments  remote  from 
each  other  can  be  quite  different.  In  the  language  of  frequency  modu¬ 
lation,  this  means  that  the  modulation  index  may  be  large.  It  is  a 
familiar  result  that  such  a  "deeply  modulated"  signal  may  have  side¬ 
bands  that  are  larger  than  the  carrier,  or  even  a  missing  carrier. 

This  is  the  phenomenon  of  spectral  line  splitting  seen  in  Figures  8,9 
and  10. 


G)0«Z03SH»I  mOC-4*-«rT>3P 


Figure  9.  Spectrum  110  and  10.5  10.5  .5 


5.  SUMMARY 


This  work  is  devoted  to  a  general  study  of  the  one -band  oscillation 
of  an  electron  under  a  high  applied  electric  field,  the  so-called  Zener 
oscillation. 

First  we  review  the  acceleration  theorem  and  the  introduction  of 
the  one-band  approximation.  The  existence  of  the  Stark  state  and  Stark 
ladder  seems  to  be  definitely  theoretically  established  in  the  litera¬ 
ture. 

We  study  general  Zener  oscillation  states  which  correspond  to 
wavepackets  of  Stark  states  or  Houston  states.  The  analysis  is  carried 
out  in  terms  of  the  expectation  values  and  uncertainties  of  both  the 
position  and  quasimomentum  operators. 

The  position  expectation  and  the  size  Ax  of  the  Zener  electron 
oscillate  in  time  with  an  amplitude  that  depends  on  the  wavepacket.  We 
discuss  in  detail  two  different  types  of  modes:  the  real  wavepacket  of 
Stark  states  with  electric-field -dependent  minimum  size;  the  real  wave- 
packet  of  Houston  states  with  arbitrary  minimum  size  and  a  resulting 
more  or  less  large  fluctuation  of  the  size  Ax. 

We  determine  that  the  Houston  state  is  the  only  state  to  satisfy 
exactly  the  equality  in  the  Heisenberg  relation:  AxAk^  >_  y  .  However, 
for  narrow  wavepackets  with  Ak^  small  the  Heisenberg  equality  is  nearly 
realized. 

We  determine  a  classical  state  corresponding  to  a  quasicoherent 
state  AxAkx  Z  y  ,  with  maximal  oscillation  amplitude  and  minimal  size 
fluctuation  of  the  Zener  electrons.  This  classical  state,  for  a  model 
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Computation  of  <^v(r,t)  lHtotalkv(r,t)> 


H.  -  =  Hrt-eex 
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APPENDIX  6.3 

Bloch  Function  Matrix  Elements 

_  _  _  2  _ 

<k' ,n|x|k,n>  and  <k' ,n'jx  |k,n> 

Differentiate  the  expression 

f  P*t  (k*  ,r)ijj  (k,r)dr  *  6  ,  <S(k-k*) 
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with  respect  to  k^: 


W-  ^n'n  =  af 

x  x 


/ei(k_k'),r  u*,(k')u  (k)dr 
j  n  n 


.  f  i(k-k’)*r  *  ,r-, .  ,r\ 

=  i  /  e  u  ,  (k  )xu  \.k)dr 

J  n  n 


i(k-kf)*t 

e 


u*,(k’) 


u  (k)dr 
n 


(C2) 


=  i<k* ,n' | x| k,n>  -  iX  ,  3(k-k?)  (C3) , 

’  1  1  n  n 

3 

defining  Xn'n  =  f  u*,  (k' ,7)  -- un(k,r)d?  . 

u.  c . 

3 

The  step  leading  to  (C3)  can  be  taken  because  u*t  has  the  crystal 

x 

periodicity  and  hence  the  last  integral  in  (C2)  is  proportional  to 
d(k-k’). 

So  we  have  the  result 

<k*,n*|x|k,n>  »  -id  ,  <5 (k-k r )  +  X  .  (k’)d(k-k’) 

1  n  n  aic  n  n 

:c 

Now  differentiate  both  (C2)  and  (C3)  with  respect  to  k^: 
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We  can  use  (C4)  and  (C5)  to  evaluate  the  troublesome  middle  integral  on 


the  rhs  of  both: 
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For  the  states  of  a  single  band,  n»n*,  we  can  take  advantage  of  the 
Hermitean  property  of  the  coordinate  operator  to  simplify  this  expres¬ 
sion.  Starting  with  the  normalization  integral 
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are  pure  imaginary.  But  since  |3un/3x|  is  real  we  can  define  i 
quantity 
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In  this  expression  the  last  term  is  imaginary,  and  it  cancels  the 
imaginary  part  of  the  third  terra  in  the  bracket. 
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For  a  narrow  wavepacket 
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<  e-*  0  using  the  Chernoff  bound  on  the  error  function.  In 
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our  calculation  —  <  •!  the  error  is  e  <  3,8  10 


DATA BAND 


90  DJHEBUO,10,10) 


100  REN  DSTA'BAND 
110  FOR  Y  s  3  TO  7 

IS  Ua  E0.95175!o,5256  >0.4920  >0,8402  >1.3774 

140  NEXT  Y 

150  FOR  Y  =  4  TO  4 

160  READ  EH5,Y,5) 

170  DATA  0.4920>0.2352  >0.4920 
180  NEXT  T 

190  EBU>5>5)  =  0,4920 
200  FOR  Y  =  2  TO  9 
210  READ  EK3>Y>4) 

220  DATA  1.2081, 0.7419,0.6448, 0.9517,1.3834, 1.4178, 1.5954, 1.5174 

230  NEXTY 

240  FOR  Y  =  3  TO  8 

250  READ  EB(4>Y>4) 

240  DATA  0.4448>0.3333>0.5254>1.0043>1,4140>1.5591 

270  NEXT  Y 

280  FOR  Y  =  4  TO  7 

290  READ  EB(5>Y>4) 

300  DATA  0,5254>0.4920>0,8402>1»3787 
302  NEXTY 
“  “  EB(6>5>4)  =  0,3402 
EK  6,6,4)  =  1.0043 
FOR  Y  =  1  TO  10 
READ  EB(2>Y>3) 

DATA  1,5422>1.0875>0,9547>1,2081>1.5525>1,4272>1.4612>1.2445>1.1244>1.0796 
...  NEXTY 
360  FOR  Y  =  2  TO  10 
370  READ  E£K3,Y,3) 

380  DATA  0,9567>0,6097>0,7419>1,1476>1.4292>1.4058>1,2594>1.1368>1.0937 
390  NEXTY 
400  F0RY  =  3T0  9 
410  READ  EfK4,Y,3) 


420  DATA  0.7419,0.6443,0.9517,1.3336,1.6178,1.5956,1.5174 


430  mV-  4 


440  FOR  Y  =  4  TO  8 
450  READ  EB(5>Y>3) 

440  DATA  0.9517,1. 0333, 1.3776>1»3322>2.1598 

470  NEXTY 

480  FOR  Y  =  5  TO  7 

490  READ  EK6,Y,3) 

500  DATA  1.3776,1.4160,1.6173 
510  NEXTY 

520  EB<  7,6,3)  =  1.6173 
530  FOR  Y  =  0  TO  10 
540  READ  EK1»Y,2) 

550  DATA  1.3548,1.2464,1.2442,1.5422,1.8326,1.7123,1.4302,1,1628,0,9516,0.8188,0.7726 

560  NEXT  Y 

570  FOR  Y  *  1  TO  10 

580  READ  EB(2,Y,2) 

590  DATA  1.2442, 0.9695, 1.0375, 1.4402, 1.58-65, 1.4180, 1.1596, 0.9330, 0.7841, 0.7316 

600  NEXTY 
410  FOR  Y  *  2  T0o10 

DATA  ^?,0sl5?0.9567» 1 ,2031 ,1.3323,1.6272,1, 4612, 1 .2645, 1 .1264,1 .0796 
NEXTY 

FOR  Y  *  3  TO  10 
READ  EB(4,Y,2) 

_  .  DATA  1.2081,1.2455, 1.3404, 1,3999,2. 0175,1. 9071, 1.7938, 1,7530 

480  NEXTY 
m  FOIUM  T029 

710  lJ404‘l.3737, 1.7734, 2, 13??, 2.6307, 2, 6853 

$0  KV*  5  TO  8 
740  READ  EH  6,  Y, 2) 

7»  DATA  1,3734,1,5591,1.3956,1.9071 
740  NEXT  Y 
770  EH  7,4,2)  =  1.5936 
780  07,7,2)  *  1.2394 
790  FOR  Y  *  0  TO  13 
MO  READ  EKOfY,l) 


BATA  0.5005,0,7479,1,3543,2,0043,1.4664,1.2703,0.9484,0.7105,0.5432,0.5015,0.4972 


FOR  Y  =  0  TO  10 
READ  EB<lrYrl> 

DATA  0.747?i0.8537rl.2444rl.759irl.7343tl.3392rl.0477r0.8107r0,4295»0.5228r0.4874 

KV*  1  TO  10 
880  REAR  EB(2rY»l) 

DATA  I,2444rl.2442rl.3422ri.3324ri.7137il.4302rl,i423r0.9514r0.8183r0.7724 
NEXT  Y 

FOR  T  *  2  TO  10 
READ  EM3rY»l) 

930  DATA  l,5422»1.5409»1.8324i2.0357>2,0010fl.7745»1.5723»1.4402»1.3717 
NEXT  Y 

FOR  r  =  3  TO  10 
READ  Eft 4>Yrl ) 

970  DATA  I,8325rl,8514r2»0667r2.4171r2.5481r2»4142r2.292?r2.2503 

NEXT  T 

FOR  Y  =  4  TO  10 
READ  EKSrYrl) 

1010  DATA  2»0447rl.8804il,9142r2.2077»2.4347j3.i454r3»i?39 
NEXT  Y 

FOR  Y  s  5  TO  9 
READ  EB(4rYil) 

1050  DATA  1.9142rl.5513rl.5174rl,7938r2.2?29 

1060  NEXT  Y 

1070  F0RY  =  4T0  8 

1060  READ  EK7rYrl) 

1090  DATA  1.5174rl,1343rl.l244 
1100  NEXT  Y 
1110  EK8r7rl)  =  1.1244 
FOR  Y  =  0  TO  10 
READ  EKOrYrO) 

•DATA  0r0.5005rl.3455rl.3974  rl.4450rl.0414r0.7353r0.4954  r0.3547r0.3018r0 

NEXT  Y 

FOR  Y  =  0  TO  10 
READ  EftlrYrO) 

1180  DATA  0,5005r0.7479rl,3548>2.0045rl.4444ri.2703r0.9484f0.7103r0.5432r0.3015r4972 

1190  NEXT  Y 

FOR  Y  =  1  TO  10 
READ  EK2rYrO) 

DATA  1 .3543r 1 ,5133r 1,9341 r 2. 1835r 1 ,3290r 1 ,5170r 1 ,2303rl.l227r 1 ,0445» 1 ,0252 

FORY  =  2  TO  10 
READ  EW3rYr0) 

DATA  1.9340rl.9745i2.2435r2,4054r2.3194r2.0343>l»9244rl,8307r 1,7992 
0  NEXT  Y 
FOR  Y  =  3  TO  10 
READ  EBUrYrO) 

DATA  2.2435r2.1542r2.2747r2.4204r2.9899r2.8444r2.7444»2.7239 
NEXT  Y 

FOR  Y  =  4  TO  10 
READ  EKSrYrO) 

DATA  2.2747rl,9534rl,9352j2.1955r2»4553r3,1939i3,4434 
NEXT  Y 
1360  FOR  Y  =  5  TO  10 
1370  READ  EftfefYrO) 

1380  DATA  1.9352rl.5341il.4850rl.7530r2.2507r2.723? 

NEXT  Y 

FDR  Y  =  4  TO  9 
READ  EB(7rYrO) 

DATA  1.4850rl.0934rl,0794rl.3930 
NEXT  Y 

EK8r7rO)  *  1,0794 
Eft  3r3r0)  =  0,7315 
1400  PRINT  CHS!  <4)f*BL0AD  CHAIN»A520* 

1410  CALL  520TILL8Z3’ 


F I LLBZ2 


390 

400 

410 

412 

414 

416 

418 

420 

422 

424 

500 

510 

520 

530 

540 

550 

560 

570 


600 

610 

620 


REN  F1LLBZ3 

PgNT  'FIL^BRJLLQUIN  ZONE* 

mills 

JFEKX»Y»Z)<  >  OTHER  190 
EHX»Y»Z)  =  EKXrZfY) 

IFEKXjY»Z)<  >  OTHER  190 
mxtbi)  =  mirhx) 

NEXT  Z 
NEXT  Y 
NEXT  X 

PRINT  'END  OF  3/48  BZ* 

FOR  X  *  0  TO  10 
FOR  Y  -  0  TO  10 
FOR  Z  =  0  TO  10 
IF  EB(X»Y»Z)  <  >  0  THEN  280 
EKX»Y>Z)  =  EB(10  -  Y»10  -  X»10  -  Z) 
NEXT  Z 
NEXT  Y 
NEXT  X 

PRINT  'END  OF  6/43  BZ* 

FOR  X  =  0  TO  10 
FOR  Y  =  0  TO  10 
FOR  Z  =  0  TO  10 
IF  EKX»Y»Z)  <  >0  THEN  380 
EBtXiYrZ)  =  EBtYfXrZ) 

NEXT  Z 
NEXT  Y 
NEXT  X 

PRINT  "END  OF  12/48=1/4  BZ  OR  1/8  CUBE 

FOR  X  =  0  TO  10 

FOR  Y  =  0  TO  10 

FOR  Z  =  0  TO  10 

IFE3(XfY.Z)<  >  0  THEN  420 

PRINT  *E8( 'XVYVZ*  )* 

NEXT  Z 
NEXT  Y 
NEXT  X 

W  =  CHR5  (4)5  REN  CTRL-D 
PRINT  MS'OPEN  CUBEBA®’ 

PRINT  WS’BELETE  CUBEBAND* 

PRINT  BIROPEN  CUBEBANB* 

PRINT  D4*"JRITE  CUBEB4MD* 

FOR  X  =  0  TO  10 
FOR  Y  =  0  TO  10 
FOR  Z  =  0  TO  10 
^MfT^EB(XiY>Z) 

leg  y 

NEXT  X 

PRINT  MS 'CLOSE  CUBEBANB* 


-100- 


D1RZENER2 


100  REM  MKCTffiTPRKRMi  . 

105  PRINT  ‘DIRM2CK4IJO' 

400  REM  PARAMETER  SECTIffit 
405  PI  -  3,141592653 
410  REM  LATTICE  PA? AM, 

415  LP  =  5,64 
420  Li  =  PI  /  LP 

503  PRINT  'INPUT  PERIODIC  VECTOR  K(«l)' 

505  INPUT  "H=*jK<0) 

504«0>  =  l«0)*2 
510  INPUT  *(<«*»*«  1) 

512  K(  1 )  s  K(  1 )  I  2 
§15  INPUT  *L=*>K(2) 

\  -  y(7 )  t  2 

520  "PRINT  'INPUT  INITIAL  VECTOR  KOItfiO.)' 

525  INPUT  'H=*»KOIO) 

““  KOI  0)  =  K010)  t  2 
INPUT  'K=*;K(K1) 

KIN^t['Ls,!k)0?2] 

KpRifiT"'Gli|,0*KR  OF  INTERPOLATION* 

INPUT  'ORDER' 50  „ 

PRINT  'GIVE  NUMBER  OF  POINTS' 

PRINT  'TO  BE  INTERPOLATED  FOR* 

PRINT  'FOURIER  SERIES  CALCULATION*  .......  M  . 

KL  =  UK<0>/  2>  t  2  f  (Ml)  /  2)  t  2  f  CKC2)  /  2)  t  2)  t  (1  / 
PRINT  'KITH  H>=*  INT  (21  *  KL) 

INPUT  'H=')H 
HOHE 

PRINT  'DATA  GIVENl*  ,  "  ,  „„ 

PRINT  *K=('K<0)  /  2VK(1)  /  2VK(2)  /  2*)* 


527 

510 

512 
515 
517 

513 
53? 
f>40 
542 
54* 
545 
5  4-5 
550 
570 
572 
575 
530 
532 
535 
595 
597 
600 
605 
610 
617 

614 

615 

616 
61? 
620 
630 
640 
650 


2) 


•KJH'KOCO)  /  2VKK1)  /  2VK0(2)  /  2*>* 
. ......  O'ORDER* 

PRINT  *M='H 
PRINT  "OK?' 

INPUT  B* 

IF  B$  =  'OK*  THEN  610 
GOTO  500 

REH  CHAIN  HITH  INTER?.  ^ 

PRINT  'IS  THE  DISK  IN  THE  DRIVE1?* 

INPUT  A* 

IF  0  =  2  THEN  640 
IF  0  =  3  THEN  620 
GOTO  533 

PRINT  CHR*  (4)}*BL0AD  CHAIM.A520' 

CALL  520'PARA3RD2*  __ 

PRINT  CHRt  (4)»*BL0AD  CHAIN  >A520* 

CALL  520'PAPA2ND2* 


REH  PARABOLIQ'JE  IHTERP.3RB 
PRINT  "PARA3PJ)  CHAINED* 

KF  FN  B<X)  =  ABSiX-5F5*(-lt 
FOR  I  =  0  TO  H  -  1 


PARA3RD2 


INT  <X  /  10))  -  10  *  IHT  (X  /  10)) 


FOR  I  =  0  TO  H  -  1 
Eli)  =  0 

xm  =  l  *  XL  /(H-  DILI 
FOR  J  =  0  TO  2 

REH  Y(J>  IS  A  POIHT  OF  INTERPOLATION 
r« J)  =  10  *  (KU)  I  I  /  <H  -  1)  F  KOU)) 

REH  KJ)  CLOSEST  POIHT 

H(J)  =  INT  (Y(J))  F  IHT  (1  /  2  F  Y(J)  -  INT  <T(J))) 

DCJ)  *  Y< J)  -  SICJ) 

NEXT  J 

HI  =  FH  B(U(0)) 

H2  =  FN  B(Hd)> 
i  i)3  =  FH  8(H!2)> 

:  E  =  E8(U1*U2*H3) 
i  IF  KO)  >  0  THEN  276 
!  A1  =  -  1 
GOTO  273 
A1  =  1 

TF  Kl)  >  0  THEN  284 
A2  =  -  1 
GOTO  284 
A2  =  1 

IF  K2)  >  0  THEH  292 
43  =  -  1 
GOTO  310 
A3  =  1 

■  M  ■  FH  B(H(0)  F  1) 

W5  =  FH  B(«d)  F  1) 
d!6  =  FNB'H'2)  FI) 

W  =  FH  8(9(0)  -  1) 

59  =  FH  BfHtl)  -1) 

51?  ■  FN  B(5K2)  -  1) 

91  =  FH  B(!nO)  F  Al) 

'.'2=  FHB55KDFA2) 

93  =  FH  8<iK2)  F  A3) 

94  =  FH  B'.SKO)  F  2  t  Al) 

',*5  =  FH  8(11(1)  F  2  I  A2) 

96=  FH  B(H'2)  F  2  *  A3) 

AD)  =  (U  I  Al  F  2)  $  (EKW4.W2*H3)  -  E)  F  (2  -  4  I  Al)  t  <EB(H7*«2»«3)  -  E)  F  E  -  EB(W*H2,H3)>  /  6  /  Al 
KO)  =  (EB5B4,H2,H3)  F  EB(W7,H2,«3))  /  2  -  E 

54 0)  =  (  -  (2  F  Al)  t  (EB(H4*H2*K3)  -  E)  F  (Al  -  2)  t  <EB(B7*H2*y3)  -  E)  -  E  F  EB(94,H2*H3)>  /  4  /  Al 

A(l)  =  ((4  t  A2  F  2)  t  (EB(H1,H5»«3)  -  E)  F  <2  -  4  I  A2)  *  <EB(H1»HS»H3)  -  E)  F  E  -  EBtHl*95*y3))  /  6  /  A2 

B(l)  =  (EBU1*H5*W3)  F  EB(W1*«S,H3))  /  2  -  E 

Gd)  =  (  -  (2  F  A2)  *  (EKK1jH5»H3)  -  E)  F  (A2  -  2)  t  <EKUfH3»H3)  -  E)  -  E  F  EB(H1*V5*H3)>  /  6  /  A2 

#2)  =  «4  I  A3  F  2)  *  (EKIUrlfifif*)  -  E)  F  (2  -  4  *  A3)  *  (EB(H1»H2>H9)  -  E)  F  E  -  EB(HliH2»W))  /  6  /  A3 


AC2)  =  ((4  $  A3  F  2)  *  (EB<W1,W2,H6)  -  E)  F  (2  -  4  *  A3)  *  (EKtfl'K 

3(2)  =  (EBCH1»U2»H4>  F  EKHl»IBfW»  /  2  -  E 

G(2)  *  (  -  (2  F  A3)  I  (EB(H1*H2*W6)  -  E)  F  (A3  -  2)  *  <EK»I*W2*H?) 

H  =  EKV1»92*H3)  -  EB(91*H2*H3)  -  EB(Hlf92*H3)  F  E 

L  =  EB(V4.92*H3)  -  EB!94*W2iH3)  -  EBU1*92*H3)  F  E 

H  =  £8(91*95*513)  -  EB(V1»H2,H3)  -  EB(W»V5*H3)  F  E 

F(2)  =  Al  *A2*(3*H-(LFN)/2) 

Hd)  =  A2  I  (L  /  2  -  H) 

K2)  =  Al  I  (H  /  2  -  H) 

1  H  =  EB(91»H2iV3)  -  EB(V1>H2*H3)  -  EB5H1jH2*93)  F  E 
L  =  EB(94*H2*V3)  -  EBU'4f«2»«3)  -  EKUlfN2»V3>  F  E 
H  =  £B(91jH2j96)  -  EB(9i*H2*K3)  -  EB(H1»H2A'6)  F  E 
F(l)  =  Al  I  A3  I  (3  f  H  -  (L  F  N)  /  2) 

:  H(3>  =  A3  t  (L  /  2  -  H) 
rt(6)  =  Al  *  CH  /  2  -  H) 

H  =  EKW*92*93)  -  EBU1*92*H3)  -  EB(W1*H2*93)  F  E 
L  =  E5( HI *92*96)  -  EB(U1*92*«3)  -  EB(Hl*y2*96>  F  E 
H  =  EBM*95*93)  -  EB(H1*H2>93)  -  EB(H1*V5*«3)  F  E 
F(0)  =  A3»A2t(3*H-(LFH)/2) 

!H(3>  =  A2  t(L/2-H) 

H(4)  *  A3  4(H  /  2  -  H) 

H4 7 )  =  EK  91*92*93)  -  EB(91*H2*H3)  -  EB(«1*92*H3)  -  EB(H1*»2*93)  F  : 
A2 

:  HI7>  =  (H(7)  -  Al  I  (H(2)  F  HC6))  -  A2  *  (Hd)  F  H(3))  -  A3  t  (H(4) 
[£(!)  =  EBU1*H2*H3) 

FOR  J  =  0  TO  2 

Ed)  =  Ed)  F  A( J)  »  IKJ)  F  8U)  *  DU)  t  2  F  GU)  t  DU)  f  3 
!  NEXT  J 


E)  -  E  F  EB(HI*IC*96))  /  4  /  A3 


EB(H1*H2*93)  F  2  »  E  -  F(0)  I  A2  «  A3  -  Fd)  I  Al  *  A3  -  F(2)  t  Al  t 
[»  -  A3  t  (HU)  F  H(5)))  /  Al  /  A2  /  A3 


504  ECI>  =  EC  I)  *  H(l)  *  BOO)  t  2  «  BC1)  +  HC2)  t  DCO)  *  DC1)  f  2  F  F(2)  *  BCO)  *  BCD 

503  E(I)  =  HI)  +  HC5)  t  BOO)  t  2  *  BO 2)  +  HC6>  *  BOO)  t  BC2)  t  2  +  F<  1 )  t  BOO)  t  BC2) 

512  HI)  =  HI)  F  FCO)  *  BCD  *  BC2)  F  HC3)  t  BCD  *  BC2)  t  2  F  HC4)  *  BCD  t  2  *  9C2) 

516  HI)  =  HI)  F  HC7)  t  BOO)  *  BCD  *  BO 2 ) 

540  PRIffT  •EC  I*  >=aE(I) 

cr.j  liryj  r 

590  PRUfT^IS  THE  DISK  IH  THE  BRIVE1?“ 

joe  njoirr 

400  PRINT  CH30  CHAIN- A520’ 

610  CALL  520,DIS?ATCH2‘ 


REH  FOURIER  CALCULATION  F0URSERIES2 

PRINT  "FO'JRSERIES  CHAINED* 
n  =  20 

PI  =  3.141592653 
A  *  KL  *  PI  /  LP 
PRINT  “COMPUTE  WI)  FRON  1 
INPUT  “I=“)EX 
PRINT  *T0* 

INPUT  "I=VFZ 
IF  il  >  0  THEN  200 
IP  FZ  >  100  THEN  160 
RIO)  =  0 

;  FOR  N  =  0  TO  H  -  2 

RIO)  =  RIO)  t  IXIN  41)-  X(N))  t  IEIN)  4  E(N  4  D)  /  2 
NEXT  N 

RIO)  =  RIO)  /  A 
EZ  =  1 

PRINT  *R(0)=*R(0) 

FOR  K  =  EZ  TO  FZ 
H  =  K  *  PI  *  2  /  A 
alX)  =  0 

:cik)  =  0 

RIK)  =  0 

FOR  H  *  0  TO  H  -  2 
L  =  N  4  1 

B  =  (E(L>-E(N))/IX(L)-X(N)) 

CD  =  I  EC N )  I  XIL)  -  EC L)  *  X(N»  /  (XIL)  -  X(N» 

C  =  CD  4  XIL)  I  D 
F  =  CD  4  XIN)  *  D 
i  CIK)  =  C(K)4Ct  SIN  (H*  XIL)) 

CIK)  =  CIK)  -  F  >  SIN  (HI  XIN)) 

CIK)  =  CIK)  4  D  /  H  1 1  COS  (H  t  XIL))  -  COS  (H  *  XIN))) 

3(K)=S!K)-CI  COS  (H*  XIL)) 

SIK)  =  SIK)  4  F  *  COSIHIXIH)) 

SCK)  =  SCK)  4  D  /  H  *  I  SIN  (H  »  XIL))  -  SIN  IH  *  XIN))) 

NEXT  N 

SCK)  =  SCK)  /  K  /  PI 
CIK)  =  CIK)  /  K  /  PI 
RIK)  -  (SIK)  t  2  4  CIK)  t  2)  t  (1  /  2) 

PRINT  “SI'K*  )=’S(K) 

.  PRINT  'CCK*  )=*C(K) 

PRINT  *R('K*)=*R(K) 

NEXT  K 

PRINT  "IS  THE  DISK  IN  THE  DRIVE  1?* 
i  INPUT  A* 

PRINT  CHR4  (4))"BL0AD  CHAIN»A520* 

CALL  520*DISPATCH2* 


-104- 


m 

102 

104 

105 
104 
103 
110 
111 
112 
140 
150 
140 
170 
130 
l?fl 
200 

210 

VC 

250 

255 

254 

260 

270 

280 

290 

300 

3lo 

320 

330 

350 

355 

340 

3<" 

370 

nr 

w  /  O 

330 

385 

390 

nr 

V.'  if 

400 

405 

410 

415 

420 

425 


RE*  DISPATCH  PROGRAM 
HOME 

PRINT  -PISPATCH  CHAINED' 

PRINT  'FOURIER  COMPONENT  COMPUTED’ 

FOR  I  =  0  TO  100 
IF  8(1)  *0  THEN  112 
PRINT  ’  R(’I*  )* 

GX  *  I 
NEXT  I 

PRINT  ’CHOOSE  THE  NEXT  PROCESS' 

PRINT  ’<1>  IS  BANBPRIHT' 

PRINT  *<2>  IS  BAM5PL0T* 

PRINT  *<3>  IS  FOURIER  SERIES  COMPUTATION' 
PRINT  “<4>  IS  FORINT' 

PRINT  ’<5>  IS  FOURPLOT* 

’<4>  IS  INPUT  NEH  K  AND  K0* 
’<7>SAVES  FOUR  IN  FOURFILE’ 

* <3> SAVES  BAND  IN  BANDFILE’ 

.....  'CHOICE!' JCH 
PRINT  "IS  THE  DISK  IN  THE  BRIVE1?' 

INPUT  AP 

THEN  350 
THEN  340 
THEN  370 
THEN  380 
THEN  390 
THEN  400 
THEN  410 
IF  CH  =  8  THEN  420 
PRINT  CHRP  <4)!'8L0AD  CHAIN»A520’ 

CALL  520’?AND?RINT2’ 

PRINT  CHRP  (4)»’BL0AD  CHAIN 
CALL  5?0'?AMBoL0T2' 

PRINT  CHRP  <4>?'BLDAD  CHAIN 
CALL  520'F0URSERIES2' 

PRINT  CHRP  (4 H'BLOAD  CHAIN 
CALL  52t}'FDimiNT2‘ 

PRINT  CHRP  U)»'BLOAD  CHAIN 
CALL  520*FGUr<PL0T2' 

PRINT  CHRP  <4>»*BL0AD  CHAIN 
CALL  520'BIRZEHER2' 

PRINT  CHRP<4)f'9L0A0CHAIN»A520’ 

CALL  520i‘FGURFILE2* 

PRINT  CHRP  (4)5*BL0AD  CHAIN»A520' 

CALL  520’BANDFILE2* 


PRINT 

PRINT 

PRINT 

INPUT 


IF  CH  =  1 
IF  CH  =  R 
IF  CH 
IP  CH  = 

IF  CH* 

IF  CH  = 

IF  CM  = 


3 

4 

5 

6 
7 


»A520‘ 

>A520‘ 

»A520* 

»A520‘ 

>A520' 


DISPATCH  2 


ZENER2F 


52  DIM  E'MO) 

54  DIM  R<  100 ) 

55  DIM  SUOO) 

53  DIM  (X 100) 

60  DIM  XdOOO) 

64  L?  =  5,64 

65  DIM  K(2) 

63  DIM  KOI  2) 

100  PRIKT  ’RETRIEVE  BANDN  H  K  L  HO  KO  LO* 

102  INPUT  ’H’iH 
10'  Ifi-3UT  *K*}K 
104  IriPOT  'L’lL 

1.05  XL  =  (H  t  2  4  K  t  2  4  L  t  2)  t  (1  /  2) 
*06  INPUT  ’HO’JHO 
107  INPUT  ’KO’fKO 

103  INPUT  ’LO’lLO 
10’  Ml)  =  2  *  K 
110  PRINT  ’GIVE  H* 

112  INPUT.  "N=’5H 

114  A*  =  'BAND*  4  STRf  (H)  4  *  * 

116  Ai  =  Al  4  STRI  (H>  4  ”  4  STRS  <K>  4  ’ 

113  Al  =  Al  4  STRI  (HO)  4  ’  ’  4  STRI  (KO)  4 
120  PRINT  ’FILE!  ’Al 

125  PRINT  ’THIS  FILE  IS’ 

130  PRINT  ’ON  WHICH  DRIVE  V 

140  INPUT  ’DRIVE  *’)DR 

141  Bl  =  STRI  (2) 

142  IF  DR  =  1  GOTO  150 

143  Bt  =  STRS  (1) 

150  D*  =  CHRS  (4)1  REM  CTRL-D 
155  Cl  =  Al  4  *»S6»D*  4  S'iRI  (DR) 

160  PRINT  DD’OPEN’CI 

170  PRINT  Dtt’READ’AI 

130  ItfUTM 

135  INPUT  I 

190  FOR  I  =  0  TO  H  -  1 

200  INPUT  K 

210  INPUT  E(I) 

220  NEXT  I 

230  PRINT  Dli’CLOSE’AI 
240  PI  =  3.141592653 

242  LI  =  PI  /  LP 

FOR  I  =  0  TO  N  -  1 
2*6  « I )  =  I  *  KL  /  (H  -  1)  *  LI 

243  NEXT  I 

250  K(Q)  =  H  »  2 

254  M2)  *  L  I  2 

256  *3(0)  *  H3  t  2 

253  KQ-;  1)  =  KO  *  2 

260  K0(2)  =  LO  t  2 

5*0  PRINT  ’IS  THE  DISK  IN  DRIVEl’ 


’4  STRI  <L)  4  * 
”4  STRI  (LO) 


510  INPUT  A* 

520  PRINT  CH3I  (4)»’BL0AD 
530  CALL  520' DISPATCH2* 


CHAIN >A520rS6fD’BI 
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